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Abstract 

Fatgraphs are multigraphs enriched with a cycHc order of the edges 
incident to a vertex. This paper presents algorithms to: (1 ) generate the 
set 3?g,„ of fatgraphs, given the genus g and the number of boundary 
cycles n; (2) compute automorphisms of any given fatgraph; (3) compute 
the homology of the fatgraph complex ^g,n- The algorithms are suitable 
for effective computer implementation. 

In particular, this allows us to compute the rational homology of the 
moduli space of Riemann surfaces with marked points. We thus compute 
the Betti numbers of Mg,„ with {2g + n) ^ 6, corroborating known results. 



1 Introduction 

This paper deals with algorithms for the enumeration of fatgraphs and their 
automorphisms, and the computation of the homology of the complex formed 
by fatgraphs of a given genus g and number of boundary components n. 

A fatgraplfl is a multigraph enriched with the assignment, at each vertex v, 
of a cyclic order of the edges incident to v. Such graphs can be "fattened" into a 
smooth punctured oriented surface, by gluing polygons along the edges in such 
a way that two adjacent edges on the polygon boundary are consecutive in the 
cyclic order at the common endpoint (see Figure [T|); an additional assignment of 
a length for each edge allows to define a conformal structure on the surface. The 
resulting Riemann surface is naturally marked, by choosing the marking points 
to be the centers of the polygons. There is thus a functorial correspondence 
between fatgraphs and marked Riemann surfaces; a fatgraph G is said to have 



^Fatgraphs have appeared independently in many different areas of mathematics: sev- 
eral equivalent definitions are known, with names such as "ribbon graphs", "cyclic graphs", 
"maps" , "dessins d'enfants" , "rotation systems" . See |21| for a comprehensive survey. 
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genus g and n boundary components if it corresponds to a punctured Riemann 
surface S G Mg_„. 

In the papers [H] and [5D], M. Kontsevicli introduced "Graph Homology" 
complexes that relate the stable homology groups of certain infinite-dimensional 
Lie algebras to various other topological objects. In particular, the "associative 
opcrad" variant of this construction results in a chain complex whose homology 
is isomorphic to the (co)homology of the moduli space of smooth Riemann 
surfaces Mg^„: the graded module underlying the complex is freely generated 
by the set 3ig,n of fatgraphs of genus g and number of boundary components n, 
endowed with the differential defined by edge contraction. 

The needed definitions and theorems about fatgraphs and their homology 
complex are briefly recalled in Section [2] the interested reader is referred to [23] 
and [H] for proofs and context. 

The bulk of this paper is concerned with finding an effectively computable 
representation of fatgraphs (see Section [J]), and presenting algorithms to: 

(1) compute automorphisms of any given fatgraph (Section!?]); 

(2) generate the set 'Rg^n of fatgraphs, given the genus g and number of bound- 
ary components n (Section [5]); 

(3) compute the homology of the fatgraph complex 3?g.„ (Section [5]). 

Note that, in contrast with other computational approaches to fatgraphs (e.g., 
[26j ) which draw on the combinatorial definition of a fatgraph, our computer 
model of fatgraphs is directly inspired by the topological definition, and the 
algorithm for enumerating elements of 3?g_„ is likewise backed by a topological 
procedure. 

Theorem l2.2| provides an effective way to compute the (co)homology of Mg^„. 
The Betti numbers of Mg „ can be computed from the knowledge of the dimen- 
sion of chain spaces Wp of the fatgraph complex and the ranks of boundary 
operators Dp] this computation can be accomplished in the following stages: 

I. Generate the basis set of W*; by definition, the basis set is the set 3lg_„ of 
oriented fatgraphs that correspond to surfaces in Mg^„. 

II. Work out the differential D : as matrices D'^p' mapping coordi- 

nates in the fatgraph basis of Wp into coordinates relative to the fatgraph 
basis of Wp-i. 

III. Compute the ranks of the matrices D^'p> . 

Stage I needs just the pair g,n as input; its output is the set of orientable 
marked fatgraphs belonging in 3ig,n. By definition, marked fatgraphs are dec- 
orated abstract fatgraphs, and the decoration is a simple combinatorial datum 
(namely, a bijection of the set of boundary cycles with the set {!,..., n}): there- 
fore, the problem can be reduced to enumerating abstract fatgraphs. With a 
recursive algorithm, one can construct trivalent Mg^„-fatgraphs from trivalent 
graphs in Mg„i.„ and Mg_i_„+i. All other graphs in Mg^„ are obtained by 
contraction of non-loop edges. 

The differential D has a simple geometrical definition: D{G) is a sum of 
graphs G", each gotten by contracting a non-loop edge of G. A simple imple- 
mentation of Stage II would just compare each contraction of a graph with p 
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edges with any graph with p—1 edges, and score a ±1 (depending on the orien- 
tation) in the corresponding entry of the matrix D^p\ However, this algorithm 
has quadratic complexity, and the large number of graphs involved makes it 
very inefficient already for Mo,5. The simple observation that contraction of 
edges is defined on the topological fatgraph underlying a marked fatgraph al- 
lows us to apply the naive algorithm to topological fatgraphs only, which cuts 
complexity down by a factor 0((n!)^). The resulting matrix is then extended to 
marked fatgraphs by the action of graph automorphism groups on the markings 
of boundary cycles. This is the variant detailed in Section [5] 

Stage III is conceptually the simplest: by elementary linear algebra, the Betti 
numbers can be computed from the rank of matrices D^p^ and the dimension of 
their domain space. The computational problem of determining the rank of a 
matrix has been extensively studied; it should be noted, however, that this step 
can actually be the most computationally burdening. 

It is worth mentioning that V. Godin [Tl] introduced a slightly different 
fatgraph complex, which computes the integral (co)homology of Mg_„; possible 
adaptation of the algorithms to this complex and an outlook on the expected 
problems is given in Section [T] 

An effective implementation (using the Python programming language |10j ) 
of the algorithms presented here is available at http : / / code . google . com/p/f atghol 
It has so far been used to compute the Betti numbers of Mg.„ for {2g + n) ^ 6. 

Results are summarized in Table [T] the values coincide with results already 
published in the literature. References are given in the closing Section [3 to- 
gether with a discussion on the implementation performance and possible future 
directions for improving and extending the algorithms. 



bo bi 62 ^3 bi 65 be 67 bs 69 &10 611 612 



Mo,3 


1 








Mo,4 


1 


2 






Mo,5 


1 


5 


6 






1 


9 


26 


24 


Mi,i 


1 










1 








Mi,3 


1 






1 


Mi,4 


1 






4 3 


M2,l 


1 




1 




M2,2 


1 




2 


1 



Table 1: Betti numbers of Mg.„ for 2g + n ^ 6. For readability, null values 
have been omitted and the corresponding entry left blank. See Section [7] for a 
discussion of these results. 



1.1 Notation 

Algorithms are listed in pseudo-code reminiscent of the Python language syntax 
(see [29j): comments in the code listings are printed in italics font. The word 
"object" is used to denote an heterogeneous composite type in commentaries to 
the code listings: for our purposes, an object is just a tuple '(ai, a2, . . ., a^Y, 
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where each of the slots Ui can be independently assigned a valueH we write 
X.ai to denote the slot of object X. Object slots are mutable, i.e., they 
can be assigned different values over the course of time. Appendix [B] gives a 
complete recap of the notation used and the properties assumed of syntax, data 
structures, and operators. 

A great deal of this paper is concerned with finding computationally-effective 
representations of topological objects; in general, we use boldface letters to 
denote the computer analog of a mathematical object. For instance, the letter 
G always denotes a fatgraph, and G its corresponding computer representation 
as a Fatgraph object. 

Finally, if ^ is a category of which X, Y are objects, we use Eilenberg's 
notation A[X, Y) for the Hom-set, instead of the more verbose Hom^(X, Y). 

2 Fat graphs and marked Riemann surfaces 

This section recaps the main definitions and properties of fatgraphs and the 
relation of the fatgraph complex to the cohomology of Mg^n- These results are 
well-known: a clear and comprehensive account is given by G. Mondello in [24| : 
the book by Lando and Zvonkin |21j provides a broad survey of the applications 
of fatgraphs and an introduction accessible to readers without a background in 
Algebraic Geometry. 

"Fatgraphs" take their name from being usually depicted as graphs with 
thin bands as edges, instead of 1-dimensional lines; they have also been called 
"ribbon graphs" in algebraic geometry literature. Here, the two names will be 
used interchangeably. 

Definition 2.1 (Geometric definition of fatgraphs). A fatgraph is a finite CW- 
complex of pure dimension 1, together with an assignment, for each vertex v, 
of a cyclic ordering of the edges incident at v. 

A morphism of fatgraphs is a cellular map f : G ^ G' such that, for each 
vertex v of G", the preimage f~^{V) of a small neighborhood y of w is a small 
neighborhood of a tree in G (i.e., f~^{V) is a contractible connected graph). 

Unless otherwise specified, we assume that all vertices of a fatgraph have 
valence at least 3. 

If G is a fatgraph, denote V{G), E{G) and L{G) the sets of vertices, unori- 
ented edges and oriented edges (equivalently called "legs" or "half-edges" ) . 

Let G be a fatgraph, and G' be the CW-complex obtained by contracting 
an edge a S E{G) to a point. If a connects two distinct vertices (i.e., a is not a 
loop) then G' inherits a fatgraph structure from G: if (a < ai < . . . < afc < a) 
and {a < a'l < ... < < a) arc the cyclic orders at endpoints of a, then the 
vertex formed by collapsing a is endowed with the cyclic order (ai < . . . < ak < 
a[ < ... < a'l^). The graph G' is said to be obtained from G by contraction 
of a. 

Contraction morphisms play a major role in manipulation of ribbon graphs. 

^This is the definition of what is usually called a "record" in Computer Science literature, 
and lacks important features of what is generally meant by "objects" in a programming 
context. However, the Python programming language only provides objects (i.e., records are 
implemented as objects with no methods), and our algorithm implementation relies on object- 
oriented programming features. We have thus decided to keep our choice of words closer to 
the actual code. 
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Figure 1: Thickening of a fatgraph into a Riemann surface. Left column: Start- 
ing fatgraph: the cychc order at the vertices is given by the orientation of the 
ambient euclidcan plane. Middle column: Thickening of the fatgraph by glu- 
ing topological disks along the boundary components. The border of a cells is 
drawn as a dotted line; each topological disk has been given a different color. 
Right column: The resulting Riemann surface with the embedded graph. Note 
that the two starting graphs would be isomorphic when considered as ordinary 
multigraphs; they are distinguished by the additional cyclic structure at the 
vertices. 

Lemma 2.1. Any morphism of fatgraphs is a composition of isomorphisms and 
contractions of non-loop edges. 

We can thus define functors V{—), E{—) and L{—) that send morphisms of 
graphs to maps of their set of vertices, (unoriented) edges, and oriented edges. 
The following combinatorial description of a fatgraph will also be needed: 

Definition 2.2 (Combinatorial definition of fatgraph). A fatgraph is a 4- 
tuple (L, (To, (Ti, CT2) comprised of a finite set L, together with bijective maps 
(To, (Ti, tT2 : L L such that: 

» (Ti is a fixed-point free involution: ~ id, and 

» (To O (T2 = (Tl . 

Lemma 2.2. Definitions \2.1\ and \2.2\ are equivalent. 

Any two of the maps (To, (Ti, (T2 determine the third, by means of the defining 
relation cto o (T2 = ci ; therefore, to give a ribbon graph it is sufficient to specify 
only two out of three maps. 

In the combinatorial description, V{G) is the set Lq of orbits of (To, E{G) is 
the set Li of orbits of (Ti, and L{G) is plainly the set L. 

There is a functorial construction to build a topological surface S{G) from 
a fatgraph G; this is usually referred to as "thickening" or "fattening" in the 
literature. 
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Lemma 2.3. There exists a functor S that associates to every fatgraph G a 
punctured Riemann surface S{G), and to every morphism f : G ^ G' a contin- 
uous map S{f) : S[G) S{G'). 

Denote by B{G) the set L/a2 of orbits of (J2: in the topological description, 
its elements are the support of 1-cycles in H^{G) that correspond under a re- 
traction to small loops around the punctures in S{G); they arc called "boundary 
cycles" of G. 

The assignment G 1-^ B{G) extends to a functor B{—)\ by Lemma [2.11 for 
any /: Gi — >■ G2 the map B[f): B{Gi) ~> B{G2) is a bijection. 

The correspondence between fatgraphs and Riemann surfaces allows us to 
give the following. 

Definition 2.3. The number of boundary cycles of a graph G is given by 
n = |i3(G)|, and is equal to the puncture number of the Riemann surface S{G). 
If S{G) has genus g and n boundary cycles, then: 

X(G) xiS{G)) = 2 - 2g - n = 2 - 2.9 - \B{G)\, (1) 

so we can define, for any fatgraph G, the genus g, as given by the relation above. 

Lemma 2.4. If G' is obtained from G by contraction of a non-loop edge, then 
G and G' share the same genus and number of boundary cycles. 

Definition 2.4. A marked fatgraph is a fatgraph G endowed with a bijection 
ly : B{G) — ;> {1, . . . , n}. The map i/ is said to be the "marking" on G. 

A morphism / : Gi ^> G2 of marked fatgraphs must preserve the marking of 
boundary cycles: 

B{Gi) B{G2) 




{1, •••,"} 

By a slight abuse of language, we shall usually omit mention of the marking 
map V and just speak of "the marked fatgraph G" . 

2.1 Moduli spaces of marked Riemann surfaces 

Fix integers g ^ 0, n > such that 2 — 2g — n < 0. Let S* be a smooth closed 
oriented surface of genus g and X = {xi, . . . ,Xn} a set of points of 5*. 

Definition 2.5. The Teichmiiller space 

7g,r, := Conf(S')/Diff°(S',n) 

is the quotient of the set of all conformal metrics on S by the set of all difFco- 
morphisms homotopic to the identity and fixing the n marked points. 

The mapping class group Tg n is the group of isotopy classes of self-difFeomorphisms 
that preserve orientation and fix marked points: 

rg,„ := Diff+(S',n)/Diff°(S',n). 
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The topological space Mg,„ := 7g,n/'^g,n is the moduli space of (smooth) 
rt-pointed algebraic curves of genus g. It parametrizes complex structures on S, 
up to diffeomorphisms that: (1 ) are homotopic to the identity mapping on 5, 
(2) preserve the orientation of S, and (3) fix the n marked points. 

The Teichmiiller space Tg_„ is an analytic space and is homeomorphic to a 
convex domain in C'^^"^"''". Since Tg.„ is an analytic variety and Tg,n acts dis- 
continuously with finite stabilizers, Mg_„ inherits a structure of analytic orbifold 
of complex dimension 3g — 3 + n. 

Since Tg.n is contractible, its equivariant (co)homology with rational coeffi- 
cients is isomorphic to the rational (co)homology of Mg,„ (see [H VII. 7. 7]). 

2.2 The fatgraph cellularization of the moduh spaces of 
marked Riemann surfaces 

An embedding of a fatgraph G is an injective continuous map t: G S, that is, 
a homeomorphism of G onto l{G) C 5*, such that the orientation on S induces 
the cyclic order at the vertices of i(G). 

Definition 2.6. An embedded fatgraph is a fatgraph G endowed with a home- 
omorphism L between S{G) and the ambient surface S, modulo the action of 
BiS°{S). 

There is an obvious action of Tg n on the set ^g,n of fatgraphs embedded 
into rt-marked Riemann surfaces of genus g. 

If confusion is likely to arise, we shall speak of abstract fatgraphs, to mean 
the topological and combinatorial objects defined in Definition 12.11 as opposed 
to embedded fatgraphs as in Definition 12.61 above. 

Definition 2.7. A metric £ on a fatgraph G is an assignment of a real positive 
number for each edge a G E{G). 

Given a metric ^ on a fatgraph G, the "thickening" construction for fatgraphs 
can be extended to endow the surface S{G) with a conformal structure depen- 
dent on Conversely, a theorem due to Jenkins and Strebel guarantees that 
a metric can be defined on each fatgraph embedded in a surface S*, depending 
uniquely on the conformal structure on S. 

Let G be a fatgraph (embedded or abstract) of genus g with n marked 
boundary components. The set A(G) = {(G, €)} of metrics on G has an obvious 
structure of topological cell; now glue these cells by stipulating that A(G') is the 
face £c( = of A(G) when G' is obtained from G by contraction of the edge a. 
The topological spaces obtained by this gluing instructions are denoted T™™*^ 
(when using embedded fatgraphs), or M™™'^ (when using abstract fatgraphs). 
The following theorem clarifies their relation to the Teichmiiller and the moduli 
space; details can be found, e.g., in [211 Section 4.1]. 

Theorem 2.1. The thickening construction induces orbifold isomorphisms: 

Call M{G) the cell in M™™*^ corresponding to an abstract fatgraph G, and 
T(G) the cell in T™™^ corresponding to an embedded fatgraph G. 
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The functorial action of n on 3lg,n induces an action on T™™*^, which 
permutes cehs T{G) by PL isomorphisms. 

Lemma 2.5. M™™'' is the quotient space of 7^°™^ by the cellular action of 
the mapping class group Tg^nj the projection homomorphism commutes with the 
isomorphisms in Theorem \2.1\ 

Lemma 2.6. The isotropy group of the cell T[G) ^ '^g°n^^ is (isomorphic 
to) the automorphism group Aut G of the abstract fatgraph G underlying G. 

The action of Tg n commutes with the face operators, so M{G) is a face of 
M{G') iff G" is obtained from G by contraction of a non-loop edge. 

2.3 Equivariant homology of Tg „ and the complex of fat- 
graphs 

Definition 2.8. An orientation of a fatgraph G is an orientation of the vector 
space QE{G), that is, the choice of an order of the edges of G, up to even 
permutations. 

Giving an orientation on G (resp. G) is the same as orienting the simplex 
A(G) (resp. T{G)). 

If G is a fatgraph with p edges, let Wg QE{G) be the 1-dimensional 

vector space generated by the wedge products ai A . . . A ap of edges of G. 
Every / G Aut G induces a map / : E{G) — >■ E{G) on the edges and thus a map 
/, : aiA. . .Attp /(ai)A. . ./\f(ap). Trivially, /*(aiA. . .Aap) = ±aiA. . .Aa^, 
depending on whether / preserves or reverses the orientation of G. 

Definition 2.9. A fatgraph G is orientable iff it has no orientation-reversing 
automorphisms. 

Form a differential complex of orientable fatgraphs as follows. 

Definition 2.10. The complex (py»,_D) of orientable fatgraphs is defined by: 

» Wp '■= @Q Wg, where G runs over orientable fatgraphs with (2g+n — I+p) 
edges; 

:» D :— ^i{—iydi, where di : Wp — !> Wp_i is given by: 

{ai A ... A Oi A ... A Up if is not a loop and 
G/ai is orientable, 
otherwise. 

Every oriented fatgraph (G, w) defines an element lug G Wg by taking the 
wedge product of edges of G in the order given by w; conversely, any ai A . . . A 
ctp G Wg defines an orientation on G by setting lo := ai < ... < ap. 

Theorem 2.2. The T g_n- equivariant homology ofTg^n with rational coefficients 
is computed by the complex of oriented fatgraphs {W^,D), i.e., there exists an 
isomorphism: 

Hl''-{7g^„,Q) = H,iW,,D). 
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Proof. The genus and number of boundary cycles will be fixed throughout, so 
for brevity, set T := Tg^^, T := Tg,„ and 7'=°™'' := T™^'\ 
By Theorem 12. 11 we have: 

==^^('J■'°'"'',Q)■ 
Recall that iJf (T™'"'', Q) can be defined as the homology of the double complex 
P* (g) Ch,('J'=°™'',Q), where is any projective resolution of Q over Q[r]. The 
spectral sequence E^^ := Hq{P^ ®r Cp) = i7,(r, Cp) abuts to (see 
VII.5 and VII.7]). 

The space T^"™*^ has, by definition, an equivariant cellularization with cells 
indexed by embedded fatgraphs of genus g with n marked boundary components. 
Let Rp be a set of representatives for the orbits of p-cells under the action of 
r. By Lemma 12. 5[ Rp is in bijective correspondence with the set of abstract 
fatgraphs having p edges, and the orientation of a cell translates directly to an 
orientation of the corresponding graph. For each geometric simplex T{G) C 
jcomb^ let Tq be its isotropy group, and let QG be the F^j-modulc consisting 
of the (Q-vector space generated by an element A on which acts by the 
orientation character: r • A = ±A depending on whether r preserves or reverses 
the orientation of the cell T{G). By Lemma l^TBl there is an isomorphism between 
Fg and Aut G; if r £ F^ reverses (resp. preserves) orientation of T{G), then the 
corresponding / £ Aut G reverses (resp. preserves) orientation on G. Therefore, 
QG and Wq are isomorphic as Aut G = F^^ modules. 

Following [SI p. 173], let us decompose (as a F-module) 

G»(T™"*,Q)= Wg; 

GeRp 

then, by Shapiro's lemma [SI III. 6. 2], we have: 

H,{T,Cp) - i/,(Fe,QG)^ HgiAutCWa). 

GeRp GeRp 

Since Aut G is finite and we take rational coefhcients, then Hq{Avit G, Wq) = 
if (J > [SI III. 10. 2]. On the other hand, if G is orientable then Aut G acts 
trivially on Wg, so: 



Ho{Aut G,Wg) 



if G has an orientation-reversing automorphism, 
Wg if G has no orientation-reversing automorphisms. 



Let R'p be the collection of all orientable fatgraphs with p edges. Substituting 
back into the spectral sequence, we see that only one column survives: 

4,0 = ^Wg = Wp, (2) 

= for aU q > 0, (3) 

In other words, E^^ reduces to the complex {EI q, d^). 

Finally, we show that the differential : Epg Ep_i q corresponds to the 
differential D: Wp — > Wp-i under the isomorphism formula ([2]); this will end 
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the proof. Indeed, we shall prove commutativity of the following diagram at the 
chain level: 



6p &P — 1 

(4) 

which implies commutativity at the homology level: 

e^e^, i/o(Aut G, Wg) — ^ ®G'eR;_, ^o(Aut G', W^c) 

ifo(r, Gp(T-'°b^ Q)) ^ Ho{T, Gp_l(T-'°^ Q)) 

=-ffo(r,9) 

whence the conclusion g = (VK*,_D). 

The vertical maps 9p, 0p_i in Q are the chain isomorphisms underlying the 
r-module decomposition Gp(T™'"'',Q) ^ ®GeR^^G- Taking the boundary 

of a cell r(G) C T™"* commutes with the F-action: 5T(t • G) = r ■ ar(G). 
Furthermore, T(G') is a cell in dT{G) iff G" is obtained from G by contraction 
of an edge; but G" is a contraction of G iff the underlying abstract fatgraphs G' 
and G stand in the same relation. Thus, the F-complexes (G*,9) and (Ty*,D) 
are isomorphic by 6**, so diagram (|4]) commutes, as was to be proved. □ 



3 Computer representation of Fatgraphs 

Although the combinatorial definition of a fatgraph (cf. Lemma lends itself 
to a computer representation as a triple of permutations — as used, e.g., in [521 
Section 2.4] — , the functions that are needed by the generation algorithms (see 
Section [S|) are rather topological in nature and thus suggest an approach more 
directly related to the concrete realization of a fatgraph. 

Definition 3.1. A Fatgraph object G is comprised of the following data: 

» A list G. vertices of Vertex objects. 

» A list G. edges of Edge objects. 

» A set G. boundary _cycles of BoundaryCycle objects. 

» An orientation G. orient. 

The exact definition of the constituents of a Fatgraph object is the subject 
of the following sections; informally, let us say that a Vertex is a cyclic list of 
edges and that an Edge is a pair of vertices and incidence positions. A precise 
statement about the correspondence of abstract fatgraphs and Fatgraph objects 
is made in Section [3751 

There is some redundancy in the data comprising a Fatgraph object: some 
of these data are inter-depcndent and cannot be specified arbitrarily. Actually, 
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Figure 2: Representation of vertices as (cyclic) lists of edge labels; vertices are 
identified by lowercase Latin letters; edge labels are depicted as roman numerals 
on a yellow square background, sitting over the edge they label. The represen- 
tation of a vertex as a list is implicitly ciliated: here we use the convention that 
the edge closest to the tail of the arrow is the ciliated one. 



all data comprising a Fatgraph object can be computed from the vertex list 
alone, as the following sections show. 

In what follows, the letters I, m and n shall denote the number of vertices, 
edges and boundary cycles: 

» I = \V{G)\ = size{G .vertices) , 

» m ~ \E{G)\ ~ size{G .edges) , 

» 71 = \B{G)\ ~ size{G .boundary.cycles) . 

For integers a and k, we use (a%fc) to denote the smallest non-negative repre- 
sentative of a mod k . 

3.1 Vertices 

We can represent a fatgraph vertex by assigning labelf[f| to all fatgraph vertices 
and mapping a vertex to the cyclically-invariant list of labels of incident edges. 
Figure [2] gives an illustration. 

Definition 3.2. A vertex together with a choice of an attached edge is called 
a ciliated vertex. The chosen edge is called the cilium. 

Definition 3.3. If w is a ciliated vertex and e is a half-edge attached to it, define 
the attachment index of e at w as the index of edge e relative to the cilium at v: 
if a is the attachment index of e at tj, then ctq takes the cilium at v onto e. 

The attachment index at a vertex is unambiguously defined for all edges 
which are not loops; the two half-edges comprising a loop have distinct attach- 
ment indices. For brevity, in the following we shall slightly abuse the definition 
and speak of the attachment index of an edge at a vertex. 

^Labels can be drawn from any finite set. In actual computer implementations, two obvious 
choices are to use the set of machine integers, or the set of Edge objects themselves (i.e., label 
each fatgraph edge with the corresponding computer representation). 
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e2 = Edge((b,0),(b,2)) 




ej = Edge((a,0),(b,l)) 




eg = Edge( (a,l), (a,2) ) 

Figure 3: Representation of fatgrapli edges. Each edge is identified with a pair of 
endpoints, where an endpoint is a vertex together with an attachment index. In 
the figure, letters a and b denote the vertices; attachment indices are computed 
by assigning index to the edge closest to the orientation arrow's tail. 

Definition 3.4. A Vertex object v = Vcrtcx{ei, . . ., Cz) is a list of the labels 
ei, . . . , of attached edges. 

Two Vertex objects are considered equal if one is equal (as a sequence) to 
the other rotated by a certain amount. 

Note that the definition of Vertex objects as plain lists corresponds to ciliated 
vertices in a fatgraph. In order to implement the cyclic behavior of fatgraph 
vertices, the requirement on equality must be imposed; equality of Vertex objects 
can be tested by an algorithm of quadratic complexity in the vertex valence. 

If D is a vertex object, let us denote numJoops{v) the number of loops 
attached to t;; it is a vertex invariant and will be used in the computation 
of fatgraph isomorphisms. Implementations of numJoops need only count the 
number of repeated edge labels in the list defining the Vertex object v. 

3.2 Edges 

Definition 3.5. An Edge object e is an unordered pair of endpoints, so defined: 
each endpoint corresponds to a 2-tuple {v,a), where v is a, vertex, and a is the 
index at which edge e appears within vertex v (the attachment index). 

It is clear how an Edge object corresponds to a fatgraph edge: a fatgraph 
edge is made of two half-edges, each of which is uniquely identified by a pair 
formed by the end vertex v and the attachment index a. In the case of loops, 
the two ends will have the form (v, a), (v, a') where a and a' are the two distinct 
attachment indices at v. 

The other_end{e, v, a) function takes as input an edge object e, a vertex t>, 
and an attachment index a and returns the endpoint of e opposite to {v,a). 

The notation Edge{{endpoints)) will be used for an Edge object comprising 
the specified endpoints. 

Figure [3] provides a graphical illustration of the representation of fatgraph 
edges as Edge objects. 
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Figure 4: Representation of fatgraph boundary cycles. Left: How the bound- 
ary cycles are represented with corners: each boundary component is identi- 
fied with the set of triplets it encloses. Therefore the boundary cycles for the 
graph are represented by the sets {(a, 0, 1), (&, 2, 0)}, {(a, 1, 2), (6, 1, 2)}, and 
{(a, 2, 0), (6, 0, 1)}. Right: Zoom around vertex a in the left picture, to show the 
three corners identified with triples (a, i, j). The indices in the triple are attach- 
ment indices, i.e., displacement relative from the ciliated edge (the one closest 
to the arrow tail) ; they bear no relation to the labels on the edges (numbers on 
the light yellow background in the left picture). 

3.2.1 Computation of the edge list 

The edge list G. edges can be computed from the list of vertices as follows. 

The total number m of edges is computed from the sum of vertex valences, 
and used to create a temporary array P of m lists (each one initially empty). 
We then incrementally turn P into a list of edge endpoints (in the form {v,a) 
where v is a vertex and a the attachment index) by just walking the list of 
vertices: P[k] is the hst [ (t)fc,0), . . ., {vk,Zk) ] where Vk (of valence z^) is the 
fc-th Vertex in G. vertices. The list G. edges is just P recast into Edge objects. 
In pseudo-code: 

m ^ (1/2) • J2veG.vertices valence(i;) 
P ^ array of m empty lists 
for V in G .vertices: 
for (a, e) in enumerate{v): 

append {v,a) to P[e] 
wrap endpoints into "Edge" objects 
G. edges <~ [ Edge{p) for p in P ] 

3.3 Boundary Cycles 

Definition 3.6. A BoundaryCycle object is a set of corners (see Figure [231) • 

A corner object C is a triple (vertex, incoming, outgoing), consisting of a 
vertex v and two indices i = C. incoming, j = C .outgoing of consecutive edges 
(in the cyclic order at v). In order to have a unique representation of any corner, 
we impose the condition that either j = i + 1, or i and j are, respectively, the 
ending and starting indices of v (regarded list). 

It is easy to convince oneself that a BoundaryCycle object corresponds to a 
boundary cycle as defined in Section [31 Indeed, if {L, 0^,01,02) is a fatgraph. 
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Figure 5: A fat graph whose two boundary cycles are comprised of ex- 
actly the same edges; however, they give rise to disjoint sets of corners: 
{(v, 2, 3), (v, 4, 5), {v, 0, 1)} versus {{v, 1, 2), {v, 3, 4), (v, 5, 0)}. 

then the boundary cycles are defined as the orbits of (T2 on the set L of half-edges; 
a (endpoint vertex, attachment index) pair uniquely identifies an half-edge and 
can thus be substituted for it. For computational efficiency reasons, we add an 
additional successor index to form the corner triple {v,i,j) so that the action 
of (72 can be computed from corner data alone, without any reference to the 
ambient fatgrapho 

Since distinct orbits are disjoint, two BoundaryCycle objects are either iden- 
tical (they comprise the same corners) or have no intersection. In particular, 
this representation based on corners distinguishes boundary cycles made of the 
same edges: for instance, the boundary cycles of the fatgraph depicted in Fig- 
ure 13.31 are represented by the disjoint set of corners {(v, 2, 3), {v, 4, 5), (v, 0, 1)} 
and{(t;,l,2),(-u,3,4),(t;,5,0)}. 

3.3.1 Computation of boundary cycles 

The procedure for computing the set of boundary cycles of a given Fatgraph 
object G is listed in Algorithm [T] The algorithm closely follows a geometrical 
procedure: starting with any corner, follow its "outgoing" edge to its other 
endpoint, and repeat until we come back to the starting corner. The list of 
corners so gathered is a boundary cycle. At each iteration, the used corners are 
cleared out of the corners list by replacing them with the special value USED, so 
that they will not be picked up again in subsequent iterations. 

Lemma 3.1. For any Fatgraph object G representing a fatgraph G, the function 
compute_boundary_cycles in Algorithm{l\has the following properties: 1) termi- 
nates in finite time, and 2) returns a list of BoundaryCycle objects that represent 
the boundary cycles of G. 

Proof. The algorithm works on a temporary array corners: as it walks along a 
boundary cycle (lines [24ll30|) . corner triples are moved from the working array 

''This is important in order to share the same corner objcets across multiple BoundaryCycle 
instances, which saves computer memory. 
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Algorithm 1 Output the set of boundary cycles of a Fatgraph object G. Input 
to the algorithm is a Fatgraph object G; the output is a list of BoundaryCycle 
objects. The special constant used marks locations in the temporary array 
corners whose contents has already been assigned to a boundary cycle. 

dcf compute_boundary_cycles{G) : 

build working array of comers 
corners <—[[(?;, i, {{i + 1) % \v\)) for i in 0, . . . , — 1 ] 



result empty list 

io ^ 
while True: 

locate the first unused corner 
for I in Iq, . . ., size{corners) — l: 
V ^ G.vertices[l] 

i ■(— first Jndex_not-Used{corners[v], io) 
if i is not None: 
exit "for" loop 
if I = size{corners) — l and i is None: 
all corners used, mission accomplished 
return result 
else: 
lo ^ I 
io <~ i 

walk the boundary cycle and record corners 
start <— {v, i) 
triples empty list 
while {v,i) ^ start or size(triplcs) — 0: 

triples. append{corners[v] [i] ) 

j ^ corners[v][i][2] 

e ^ v[j] 
mark location as "used" 

corners[v][i] used 

{v,i) other_end{e, v, j) 

b BoundaryCycle{triples) 
result. append{b) 

def first _index_not_used{L): 

for index, item in enumerate(i): 
if item is not used: 
return index 
return None 



for V in G. vertices ] 
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to the triples list and replaced with the constant USED; when we're back to the 
starting corner, a BoundaryCycle object b is constructed from the triples list 
and appended to the result. 

The corners variable is a list, the n-th item of which is (again) a list holding 
the corners around the n-th vertex (i.e., G.vertices[n\), in the order they are 
encountered when winding around the vertex. By construction, corners[v] [i] 
has the the form {v,i,j) where j is the index following i in the cyclic order, 
i.e., {v,i,j) represents the corner formed by the "incoming" i-th edge and the 
"outgoing" j-th edge. 

The starting corner for each walk along a boundary cycle is determined 
by scanning the corners list (lines [T0Hl4|) : loop over all indexes v, i in the 
corners list, and quit looping as soon as corners[t;] [i] is not used (line [T5)) . If all 
locations in the corners list are USED, then the all corners have been assigned 
to a boundary cycle and we can return the result list to the caller. □ 

3.4 Orientation 

According to Definition 12.81 orientation is given by a total order of the edges 
(which directly translates into an orientation of the associated orbifold cell) . 

Definition 3.7. The orientation G. orient is a list that associates each edge 
with its position according to the order given by the orientation. Two such lists 
are equivalent if they differ by an even permutation. 

If ei and 62 arc edges in a Fatgraph object G, then ei precedes 62 iff 
G.orient[ei] < G.orient[e2]; this links the fatgraph orientation from Defini- 
tion with the one above. 

If a Fatgraph object is derived from another Fatgraph instance (e.g., when 
an edge is contracted) , the resulting graph must derive its orientation from the 
"parent" graph, if we want the edge contraction to correspond to taking cell 
boundary in the orbicomplex M™"*. 

When no orientation is given, the trivial one is (arbitrarily) chosen: edges 
are ordered in the way they are listed in the G. edges list, i.e., G.orient[e] is the 
position at which e appears in G. edges. 

According to Dcfinition l2.9l a fatgraph is oricntablc iff it has no orientation- 
reversing automorphism. The author knows of no practical way to ascertain if a 
fatgraph is orientable other than enumerating all automorphisms and checking 
if any one of them reverses orientation: 

def is_oriented{G): 

for a in automorphisms{G): 
if is-orientation^eversing{a) : 
return False 

no orientation— reversing automorphism found, G is orientable 
return True 

3.5 A category of Fatgraph objects 
3.5.1 Isomorphisms of Fatgraph objects 

In this section, we shall only give the definition of Fatgraph isomorphisms and 
prove the basic properties; the algorithmic generation and treatment oi Fatgraph 
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isomorphisms is postponed to Sectional 

Definition 3.8. An isomorphism of Fatgraph objects Gi and G2 is a triple 
/ ={pv, rot, pc) where: 

» pv is a permutation of the vertices: vertex Vi of Gi is sent to vertex pv[v] 
of G2, and rotated by rot[u] places leftwards; 

» pe is a permutation of the edge labels: edge e in Gi is mapped to edge 
pe[e] in G2. 

The adjacency relation must be preserved by isomorphism triples: if Vi and V2 
are endpoint vertices of the edge e, then pv[t;i] and pv[7;2] must be the endpoint 
vertices of edge pe[e] in G2. 

Since a vertex in a Fatgraph instance is essentially the list of labels of edges 
attached to that vertex, we can dually state the compatibility condition above 
as requiring that, for any vertex v in Gi. vertices and any valid index j of an 
edge of t), we have: 

G2. vertices[pv[t;]] [j+rot[t;]] = pc[Gi.vcrticcs[v][j]] (5) 

The above formula ([5]) makes the parallel between Fatgraph object isomorphisms 
and fatgraph maps (in the sense of Definition 12. ip explicit. 

Lemma 3.2. Let Gi, G2 he fatgraphs, represented respectively by Gi and G2. 
Every isomorphism of fatgraphs f : Gi —> G2 lifts to a corresponding isomor- 
phism f — (pv, rot, pe) on the computer representations. Conversely, every 
triple (pv, rot, pc) representing an isomorphism between the Fatgraph instances 
induces a (possibly trivial) fatgraph isomorphism between Gi and G2. 

Proof. Every isomorphism / : Gi — > G2 naturally induces bijective maps fv '■ V{Gi) — > 
V{G2) and fs '■ E{Gi) — > E{G2) on vertices and edges. Given a cilium on every 
vertex, / additionally determines, for each vertex v 6 V{G), the displacement 
fmt{v) of the image of the cilium of v relative to the cilium of fvi^)- Similarly, 
fE determines a bijective mapping of edge labels, and is completely determined 
by it. This is exactly the data collected in the triple (pv, rots, pc), and the 
compatibility condition ([5]) holds by construction. 

Conversely, assume we are given a triple (pv, rots, pe), representing an iso- 
morphism of Fatgraph instances. We can construct maps fv, fs as follows: fv 
sends a vertex v G Gi to the vertex corresponding to pv[t>]; fE maps the cilium 
of V to the edge attached to pv[v] at rot[v] positions away from the cilium; the 
compatibility condition ([5]) guarantees that fE is globally well-defined. □ 

Lemma 3.3. Let Gi , G2 be Fatgraph objects, and rj a bijective map between 
Gi. edges and G2.edges that preserves the incidence relation. Then there is a 
unique Fatgraph isomorphism f that extends rj (in the sense that f.pe = rj). 

Proof. Start constructing the Fatgraph morphism / by setting /.pe = rj. If 
ei, . . . , Bzf. are the edges incident to G Gi.vertices, then there is generally 
one and only one endpoint common to edges r]{ek); define f.pv[vk] = f^. 
There is only one case in which this is not true, namely, if all edges share the 
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Algorithm 2 Construct a new Fatgraph object G' obtained by contracting the 
edge e in G. Tlie renumbering function s is tlie identity on numbers in the range 
0, . . . , e — 1, and shifts numbers in range e + 1, . . . , m down by 1. Function 
rotated{L,p) returns a copy of list L shifted leftwards by p places. 

def contract{G , e): 

let (vi, ai), {v2, 02) be the endpoints of e 
V' <— [ Vcitcx{x for X iivv \i X ^ e) 

for V in G. vertices ii v Vi and v ^ V2 ] 
append the fused vertex at end of list V 
v' <— Vertcx(rotated(i»i , a{) + rotatcd{v2, 02)) 
V .append{v') 

oj' ^ [ s{G' .onent[x]) for x in G' .edges if .t 7^ e ] 
return Fatgraph{vertices <— V'; orient ^ ui') 



same two endpointslf] in this case, however, there is still only one choice of / 
.pv[vk\ such that the cyclic order of edges at the source vertex matches the cyclic 
order of edges at the target vertex. Finally, choose f.rot[vk] as the displacement 
between the cilium at v'f. and the image of the cilium of Vk ■ 

It is easy to check that eq. (O holds, so / is a well-defined isomorphism. □ 

3.5.2 Contraction morphisms 

Recall from the definition in Section [5] that contraction produces a "child" fat- 
graph from a "parent" fatgraph and a chosen regular (i.e., non-looping) edge. 

The Fatgraph.contract method (see Algorithm [2]) thus needs only take as 
input the "parent" graph G and the edge e to contract, and produces as output 
the "child" fatgraph G'. The contraction algorithm proceeds in the following 
way: 

» The two end vertices of the edge e are fused into one: the list G' .vertices 
is built by copying the list G .vertices, removing the two endpoints of e, 
and adding the new vertex (resulting from the collapse of e) at the end. 

» Deletion of an edge also affects the orientation: the orientation G' .orient 
on the "child" fatgraph keeps the edges in the same order as they are in the 
parent fatgraph. However, since G' .orient must be a permutation of the 
edge indices, we need to renumber the edges and shift the higher-numbered 
edges down one place. 

» The "child" graph G' is constructed from the list G' .vertices and the de- 
rived orientation G'. orient; the list of "new" edges is constructed according 
to the procedure given in Section 13.2.11 

Listing [5] summarizes the algorithm applied. 

The vertex resulting from the contraction of e is formed as follows. Assume 
Vi and V2 are the endpoint vertices of the contracted edge. Now fuse endpoints 
of the contracted edge: 

'''So there are only two vertices in total, and the corresponding fatgraph belongs in 3?o,m- 
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(1 ) Rotate the lists Vi, V2 so that the given edge e appears last in Vi and first 
in V2. 



(2) Form the new vertex v by concatenating the two rotated Hsts (after ex- 
punging vertices Vi and 1)2). 

Note that this changes the attachment indices of all edges incident to Vi and 
V2, therefore the edge list of G' needs to be recomputed from the vertex list. 

The "child" fatgraph G' inherits an orientation from the "parent" fatgraph, 
which might differ from its default orientation. Let ai, . . . , a/i, . . . , a„i be the 
edges of the parent fatgraph G, with e = being contracted to create the 
"child" graph G' . If ak(i) < ctk{2) < ■ ■ ■ < ak{m) is the ordering on E{G) 
that induces the orientation on G and h ~ k(j), then ak{i) < . . . < ak{j-i) < 
ak{j+i) < . . . < ak{m} descends to a total order on the edges of G' and induces 
the correct orientationl^ 

Orientation is represented in a Fatgraph object as a list, mapping edge labels 
to a position in the total order; using the notation above, the orientation of G is 
given by w := k~^. The orientation on G" is then given by w' defined as follows: 



This corresponds exactly to the assignment in Algorithm [2] 
The above discussion can be summarized in the following. 

Lemma 3.4. IfG andC represent fatgraphs G andG', andG ~ contract{G' , e), 
then G is obtained from G' by contraction of the edge e represented by e. 

The contract_boundary_cycle function. The boundary cycles of the "child" 
Fatgraph object G' can also be computed from those of G. The implementation 
(see Listing [T]) is quite straightforward: we copy the given list of corners and 
alter those who refer to the two vertices that have been merged in the process 
of contracting the specified edge. 

Let Vi and V2 be the end vertices of the edge to be contracted, and ai, 02 be 
the corresponding attachment indices. Let zi and Z2 be the valences of vertices 
Vi, V2- We build the list of corners of the boundary cycle in the "child" graph 
incrementally: the b' lists starts empty (line|6]), and is then added corners as 
we run over them in the loop between lines [7] and 1261 

There are four distinct corners that are bounded by the edge e to be con- 
tracted; denote them by Ci, C2, C3, C4. These map onto two distinct corners 
C, C' after contraction. Assume that Ci and C2 map to C: then Ci and C2 
lie "on the same side" of the contracted edge, i.e., any boundary cycle that in- 
cludes Ci will include also C2 and viceversa. (Sec Figure [H] for an illustration.) 

^That is to say, the orientation that corresponds to the orientation induced on the cell 
A(G') as a face of A(G). 
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if cj(i) < h, 
if uj{i) > h. 



Alternatively we can write: 
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Algorithm 3 Return a new BoundaryCycle instance, image of b under the 
topological map that contracts the edge with index e. 



def contract_boundary_cycle{G, b, e): 

let {vi, oi), {v2, be the endpoints of e 
zi ■(— valence{vi) 
Z2 ■(— valence{v2) 
"child" boundary cycle b' starts off as an empty list 

b'^U ^ 

for corner in b: 
if co™er[0] = Vi: 
if oi = corner. incoming: 

continue with next corner 
else: 

ii ■<— (corner. incoming — ai — 1) % zi 

«2 ^ (corner. outgoing — oi — 1) % zi 

append corner (vi, ii, 12) to b' 
elif corncr[0] = V2: 

if 02 = corner, incoming: 

continue with next corner 
if 02 = corner. outgoing: 

append (vi, zi + Z2 — 3, 0) to b' 
else: 

11 ■(— zi — 1 + ((corner.incoming — 02 — 1) % Z2) 

12 ■(— Zi — 1 + ((corner. outgoing — a2 — 1) % Z2) 
append (wi,«i,i2) to b' 

else: 

fceep comer unchanged 
append corner to b' 
return Boundary Cycle(b') 
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Figure 6: How corners are modified by edge contraction. Left: Four distinct 
corners are formed at the endpoints Vi, V2 of edge e, which is to be contracted: 
Ci = (-i;i,0,2), C2 = (-«2,0,1), C3 = (t;2,l,2), and C4 = (-^1,0,1). Edges 
are shown thickened, and (potentiaUy) distinct boundary cycles are drawn in 
different colors. Right: After contraction of e, corners Ci and C2 are fused into 
C = {v, 0, 1), and C3, C4 are fused into C = {v, 2, 3). 

Since they both map to the same corner C in the "child" graph, we only need to 
keep one: we choose to keep (and transform) the corner that has the contracted 
edge at the second index flines iQl llOp : similarly for C3 and C4 in mapping to 
C (lines HU-ini). 

Recall that, when contracting an edge with endpoints Vi and V2, the new 
vertex is formed by concatenating two series of edges: (1 ) edges attached to the 
former Vi, starting with the successor (in the cyclic order) of the contracted 
edge; (2) edges attached to the former V2, starting with the successor of the 
contracted edge. Therefore: 

(1) The image of a corner rooted in vertex Vi will have its attachment indices 
rotated leftwards by ai + 1 positions: the successor of the contracted 
edge has now attachment index (lines [T^ [T^ . Note that the highest 
attachment index belonging into this group is zi — 2: position zi — 1 
would correspond to the contracted edge. 

(2) The image of a corner rooted in vertex V2 has its attachment indices 
rotated leftwards by 02 + 1 positions, and shifted up by zi — 1 (lines \2T\ - 
[^ . As a special case, when the contracted edge is in second position we 
need to map the corner to the corner having attachment index in second 
position fline fTO)) . 

Any other corner is copied with no alterations fline l26p . 
3.5.3 The category of Fatgraph objects 

We can now formally define a category of Fatgraph objects and their morphisms. 

Definition 3.9. "X^ is the category whose objects arc Fatgraph objects, and 
whose morphisms are compositions of Fatgraph isomorphisms (as defined in 
Section m and edge contraction maps. 

More precisely, if G and G' are isomorphic Fatgraph objects, then the mor- 
phism set 3i^(G, G') is defined as the set of Fatgraph isomorphisms in the sense 
of Sectional otherwise, let m and m' be the number of edges of G, G' , and set 
k := m — to': each element in 3?'^(G, G') has the form a' o (tti o • • • o tt^) o a 
where a, a' arc automorphisms of G, G' and tti, iTk are non-loop edge 
contractions. 
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Figure 7: Construction of a fatgraph out of a set of Vertex instances: half-edges 
tagged with the same (numeric) label are joined together to form an edge. 

Theorem 3.1. There exists a functor K from the category of Fatgraph 
objects to the category 31 of abstract fatgraphs, which is surjective and full. 

Proof. Given a Fatgraph G, its constituent Vertex objects determine cyclic 
sequences vq = (ej], e?, . . . , e°J, . . . , d; = (e[,, . . . , e^,J, such that 

{eg, . . . , e°g, eJ, . . . , e!,^^^^ , Cq, . . . , e'^.^ } = {0, . . . , m — 1}. 

Fix a starting element for each of the cyclic sequences vq, . . . , vi. Then set: 

L := {{e,i,v) : v = vj G {vq, . . .,vi}, e ^ ej E v}, 

and define maps co, ci, (T2 : L ^ L as follows: 

» ctq sends {e,i,Vj) to {e' ,i' ,Vj) where i' = [i + ^)%Zj and e' = e^, is the 
successor of e in the cyclic order at Vj ; 

» (Ti maps (e, i, v) to the unique other triplet (e', i', v') G L such that e = e'; 

» finally, (72 is determined by the constraint o a2 = ci . 

Then K{G) — (i, ctq, cri, 0-2) is a fatgraph. Figure 15.5.31 provides a graphical 
illustration of the way a Fatgraph object is constructed out of such combinatorial 
data. 

Now let G be an abstract fatgraph; assuming G has m edges, assign to each 
edge a "label", i.e., pick a bijective map e : E{G) — > E, where E is an arbitrary 
finite set. Each vertex v G V{G) is thus decorated with a cyclic sequence of 
edge labels; the set of which determines a Fatgraph object G; it is clear that 
G = K{G). 

This proves that K is surjective; since every fatgraph morphism can be writ- 
ten as a composition of isomorphisms and edge contractions (Lemma 12. ip , it is 
also full. It is clear that every edge contraction is the image of an edge contrac- 
tion in the corresponding Fatgraph objects, and the assertion for isomorphisms 
follows as a corollary of Lemma 13.21 □ 

Definition 3.10. If G = K{G) then we say that the Fatgraph object G repre- 
sents the abstract fatgraph G. 
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It is clear from the construction above that there is a considerable amount 
of arbitrary choices to be made in constructing a representative Fatgraph] there 
are thus many representatives for the same fatgraph, and different choices lead 
to equivalent Fatgraph objects. 

Lemma 3.5. Two distinct Fatgraph objects representing the same abstract fat- 
graph are isomorphic. 

Proof. Assume Gi and G2 both represent the same abstract fatgraph G = 
K{Gi) ~ K{G2). Let 771, ri2 be the maps that send Edge objects in Gi, G2 
to the corresponding edges in G; then r/ = 7]^^ o rj2 maps edges of Gi into 
edges of G2 and respects the incidence relation, therefore it is the edge part of 
a Fatgraph isomorphism by Lemma 13.31 □ 

Theorem 3.2. The categories 31"^ and Ji are equivalent. 

Proof. The functor K is surjective and full by Theorem 13. II that it is also faith- 
ful follows from the following argument. Any fatgraph morphism is a compo- 
sition of edge contractions and isomorphisms. Any isomorphism determines, in 
particular, a map on the set of edges, and there is one and only one Fatgraph iso- 
morphism induced by this map (Lemma 13. 3p . Any edge contraction is uniquely 
determined by the contracted edge: if f : Gi ^ G2 is the morphism contracting 
edge e and Gi = K{Gi), then /, contraction of the Edge object e representing 
e, is the sole morphism of Gi into G2 that maps onto /. □ 

4 Fatgraphs isomorphism and equality testing 

The isomorphism problem on computer representations of fatgraphs consists 
in finding out when two distinct Fatgraph instances represent isomorphic fat- 
graphs (in the sense of Definition 12. 1[) or possibly the same fatgraph. Indeed, 
the procedure for associating a Fatgraph instance to an abstract fatgraph (see 
Theorem 13. ip involves labeling all edges, choosing a starting edge (cilium) on 
each vertex and enumerating all vertices in a certain order; for each choice, we 
get a different Fatgraph instance representing the same (abstract) fatgraph. 

The general isomorphism problem for (ordinary) graphs is a well-known 
difficult problem. However, the situation is much simpler for fatgraphs, because 
of the following property. 

Lemma 4.1 (Rigidity Property). LetGi, G2 be connected fatgraphs, and f : Gi — >■ 
G2 an isomorphism. For any vertex v G V(Gi), and any edge x incident to v, 
f is uniquely determined (up to homotopies fixing the vertices of Gi) by its 
restriction to v and x. 

In particular, an isomorphism of graphs with ciliated vertices is completely 
determined once the image w = f{v) of a vertex v is known, together with the 
displacement (relative to the cyclic order at w) of the image of the cilium of v 
relative to the cilium of the image vertex w. 

Proof. Consider / as a CW-complex morphism: / = (/o,/i) where is a 
continuous map on the set of z-dimensional cells. 

Let f/ be a small open neighborhood of S V{Gi). Given f\ij, incrementally 
construct a CW-morphism /' : Gi — > G2 as follows. Each edge x' incident to v 
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can be expressed as x' = a^x for some ^ a <Yalence(v). Let w = /(w) and 
y = /(a;), and define: 

■■= y, 

/o(w) := w, 

f[{x'):=a^y = a^f[{x) \ix' = a^x, 

foK)--=<, 

where: 

> is the endpoint of x' ~ (TqX "opposite" to v, 

> is the endpoint of y' ~ a^y "opposite" to w. 

Then /' extends / on an open set U' 3 U , which contains the subgraph formed 
by all edges attached to v and v' . In addition: 

> fi{x') = ,f{x') up to a homotopy fixing the endpoints since / commutes 
with (Tq, 

» fli{va) = f{va) since / preserves adjacency. 

By repeating the same construction about the vertices v'^ and w'^, one can 
extend /' to a CW-morphism that agrees with / on an open set U" 13 U' . 

Recursively, by connectedness, we can thus extend /' to agree with / (up to 
homotopy) over all of Gi . □ 

4.1 Enumeration of Fatgraph isomorphisms 

The stage is now set for presenting the algorithm to enumerate the isomorphisms 
between two given Fatgraph objects. Pseudo-code is listed in Algorithm |4l as 
this procedure is quite complex, a number of auxiliary functions have been used, 
whose purpose is explained in Section 14.1.11 Function isomorphisms, given 
two Fatgraph objects Gi and G2, returns a list of triples (pv, rot, pc), each of 
which determines an isomorphism. If there is no isomorphism connecting the 
two graphs, then the empty list [ ] is returned. 

By the rigidity lemma 231 any fatgraph isomorphism is uniquely determined 
by the mapping of a small neighborhood of any vertex. The overall strategy of 
the algorithm is thus to pick a pair of "compatible" vertices and try to extend 
the map as in the proof of of lemma 14.11 

We wish to stress the difference with isomorphism of ordinary graphs: since 
an isomorphism / is uniquely determined by any pair of corresponding vertices, 
the initial choice of candidates v, f{v) cither yields an isomorphism or it does 
not: there is no backtracking involved. 

Since the isomorphism computation is implemented as an exhaustive search, 
it is worth doing a few simple checks to rule out cases of non-isomorphic graphs 
(lines HBl). One has to weigh the time taken to compute a graph invariant 
versus the potential speedup obtained by not running the full scan of the search 
space; experiments run using the Python code show that the following simple 
invariants already provide some good speedup: 

» the number of vertices, edges, boundary cycles; 
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Algorithm 4 Enumerate isomorphisms between two Fatgraph objects Gi 
and G2'- output of the algorithm is a list of triples (pv, rot, pe). If there is 
no isomorphism connecting the two input fatgraphs, the empty list is returned. 



def isomorphisms{Gi, G2): 

immediately rule out easy cases of no isomorphisms 
if graphs invariants differ: 

return [ ] 
result ^ [ ] 

vsl ^ valcncc_spcctrum{Gi) 

vs2 ^ valence-spcctrum{G2) 

{valence, vertices) •(— startmg_vertices{G2) 

Vi vsl [vaience] [0] 

for V2 in compatible-vertices{vi, vertices): 
for rot in 0, . . . ,valence: 
Initialize pv, rots, pe as empty maps 

pv[vi] ^ V2 

rots[vi] ^ rot 

extend_map{pc, Vi, rotated{v2, rot)) 
if extension failed: 

continue with next rot 

breadth— first search to extend the mapping over corresponding vertices 
nexts i- ncighbors{pv, pe, Gi, Vi, G2, V2) 
while size{pv) < Gi. num. vertices: 
neighborhood -(r- [ ] 
for {v[, V2, r) in nexts: 

{pv, rots, pe) <— extendJso{pv, rots, pe, Gi, v[, r, G2, v'2) 
if cannot extend: 

exit "while" loop and continue with next rot 
append neighbors{pv, pe, Gi, v'^, G2, V2) to neighborhood 
nexts neighborhood 
isomorphism found, record it 
result. append{{pv, rots, pe)) 
return result 
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» the total number of loops; 
» the set of valences; 

» the number of vertices of every given valence. 

Since an isomorphism is uniquely determined by its restriction to any vertex, 
one can restrict to considering just pairs of the form (vi, V2) where Vi is a chosen 
vertex in Gi. Then the algorithm tries all possible ways (rotations) of mapping 
Vi into a compatible vertex V2 in G2. The body of the inner loop (line [TT] 
onwards) mimics the construction in the proof of Lemma |4. II 

The starting vertex Vi should be selected so to minimize the number of 
mapping attempts performed; this is currently done by minimizing the product 
of valence and number of vertices of that valence on G2 (line[S]), and then picking 
a vertex of the chosen valence in Gi as Vi (line[9])l3| 

First, given the target vertex V2 and a rotation rot, a new triple {pv^rots.pe) 
is created; pv is set to represent the initial mapping of Vi onto 1^2, rotated 
leftwards by rot positions, and pe maps edges of Vi into corresponding edges 
of the rotated V2- If this mapping is not possible (e.g., Vi has a loop and V2 
does not, or not in a corresponding position), then the attempt is aborted and 
execution continues from line [TT] with the next candidate rot. 

The mapping defined by {pv,rots,pe) is then extended to neighbors of the 
vertices already inserted. This entails a breadth-first searcfH over pairs of cor- 
responding vertices, starting from Vi and V2. Note that, in this extension step, 
not only the source and target vertices, but also the rotation to be applied is 
uniquely determined: chosen a vertex v'l connected to Vi by an edge e, there 
is a unique rotation r on V2 such that pv[e\ has the same attachment index 
to t>2 that e has to Vi. If. at any stage, the extension of the current triple 
(pv, rots, pc) fails, it is discarded and execution continues from line [TT] with the 
next value of rot. 

When the loop started at line [10] is over, execution reaches the end of the 
isomorphisms function, and returns the (possibly empty) list of isomorphisms 
to the caller. 

Theorem 4.1. Given Fatgraph objects Gi, G2, Junction isomorphisms returns 
all Fatgraph isomorphisms from Gi to G2. 

Proof. Given an isomorphism / : Gi G2, restrict / to the starting vertex Vi: 
then / will be output when Algorithm [4] examines the pair Vi, /(fi); since 
Algorithm [3] performs an exhaustive search, / will not be missed. 

Conversely, since equation ([5]) holds by construction for all the mappings 
returned by isomorphisms, then each returned triple / =(pv, rots, pc) is an 
isomorphism. □ 

4.1.1 Auxiliary functions 

Here is a brief description of the auxiliary functions used in the listing of Algo- 
rithm [4] and [5] Apart from the neighbors function, they are all straightforward 

'^Thc checks already performed ensure that Gi and G2 have the same "valence spectrum", 
so Gi has at least one vertex of the chosen valence. 

®The variables nexts and neighborhood play the role of the FIFO list in the usual for- 
mulation of breadth-first search: vertices are added to neighborhood during a loop, and the 
resulting list is then orderly browsed (as nexts) in the next iteration. 



26 



Algorithm 5 Enumerate the candidate extensions of the given pv and pe in 
the neighborhood of input vertices Vi and V2. 

def neighbors{pv, pc, Gi, Vi, G2, V2): 

result <— [ ] 

for each non-loop edge e attached to Vi: 
let {v'i,ai) be the endpoint of e distinct from Vi 
if v'l already in pv domain: 

continue with next e 
let (^2,02) be the endpoint of e'=pe[e] distinct from V2 
if v'2 already in pv image: 

continue with next e 
result. append{{v[, v'2, ai — 02)) 
return result 



to implement, so only a short specification of the behavior is given, with no ac- 
companying pseudo-code. 

The neighbors function. 

Definition 4.1. Define a candidate extension as a triplet {v[, t>2, r), where: 

» is a vertex in Gi, connected to Vi by an edge e; 

» t)2 is a vertex in G2, connected to V2 by edge e' = pc[e]; 

» r is the rotation to be applied to V2 so that edge e and e' have the same 
attachment index, i.e., they are incident at corresponding positions in v[ 
and v'2. 

Function neighbors lists candidate extensions that extend map pv in the 
neighborhood of given input vertices Vi (in the domain fatgraph Gi) and V2 (in 
the image fatgraph G2). It outputs a list of triplets {v'l, v'2,r), each representing 
a candidate extension. 

A sketch of this routine is given in Algorithm [5] Two points are worth of 
notice: 

(1) By the time neighbors is called (at lines [T51 and in Algorithm 2]) , the 
map pe has already been extended over all edges incident to Vi , so we can 
safely set e' =pe[e] in neighbors. 

(2) Algorithm 2] only uses neighbors with the purpose of extending pv and 
pe, so neighbors ignores vertices that are already in the domain or image 
of pv. 

The valence_spectrum function. The auxiliary function valence_spectrum, 

given a Fatgraph instance G, returns a mapping that associates to each valence 
z the list Vz of vertices of G with valence z. 
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The starting_vertices function. For each pair {z,Vz) in the valence spec- 
trum, define its intensity as the product z ■ \Vz\ (valence times the number of 
vertices with that valence). The function starting_vertices takes as input a 
Fatgraph object G and returns the pair {z,Vz) from the valence spectrum that 
minimizes intensity. In case of ties, the pair with the largest z is chosen. 

The compatible and compatible_vertices functions. Function compat- 
ible takes a pair of vertices Vi and V2 as input, and returns boolean True iff 
Vi and V2 have the same invariants. (This is used as a short-cut test to aban- 
don a candidate mapping before trying a full adjacency list extension, which is 
computationally more expensive.) The sample code uses valence and number of 
loops as invariants. 

The function compatible_vertices takes a vertex v and a list of vertices 
L, and returns the list of vertices in L that are compatible with v (i.e., those 
which V could be mapped to). 

The extend_map and extend_iso functions. The extend_map function 
takes as input a mapping pe and a pair of ciliated vertices Vi and V2, and alters 
pe to map edges of Vi to corresponding edges of V2: the cilium to the cilium, 
and so on: pe[crQ (e)] ~ ctq (pc[e]). If this extension is not possible, an error is 
signaled to the caller. 

The extend_iso function is passed a (pv, rots, pc) triplet, a vertex v[ of 
Gi, a vertex v'2 of G2 and a rotation r; it alters the given {pv,rots,pc) triple by 
adding a mapping of the vertex v'^ into vertex v'2 (and rotating the target vertex 
by r places rightwards). If the extension is successful, it returns the extended 
map (pv, rot, pe); otherwise, signals an error. 

4.2 Operations with Fatgraph Isomorphisms 

Compare pull-back orientation. The compare_orientations function takes 
an isomorphism triple (pv, rots, pc) and a pair of Fatgraph objects Gi , G2, and 
returns or —1 depending on whether the orientations of the target Fatgraph 
pulls back to the orientation of the source Fatgraph via the given isomorphism. 

Recall that for a Fatgraph object G, the orientation is represented by a map- 
ping G. orient that associates an edge e with its position in the wedge product 
that represents the orientation; therefore, the pull-back orientation according 
to an isomorphism (pv, rots, pe) from G to G' is simply given by the map 
e I— > G' .oricnt[pc[e]\. Thus, the comparison is done by constructing the permu- 
tation that maps G.oricnt[e] to G' .oricnt[pc[e]] and taking its sign (which has 
linear complexity with respect to the number of edges) . 

The is_orientation_reversing function. Determining whether an automor- 
phism reverses orientation is crucial for knowing which fatgraphs are orientable. 
Function is_orientation_reversing takes a Fatgraph object and an isomor- 
phism triple (pv, rots, pc) as input, and returns boolean True iff the isomor- 
phism reverses orientation. This amounts to checking whether the given orienta- 
tion and that of the pull-back one agree, which can be done with the comparison 
method discussed above. 
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Algorithm 6 Function MgnGraphs returns all connected fatgraphs having 
prescribed genus g and number of boundary cycles n. Actual output of the 
function is a list R, whose fc-th element R[k] is itself a list of graphs in 3^g.„ 
with m ~ k edges. 

def MgnGraphs{g ,n): 

m •<— 4^ + 2n — 5 maximum number of edges 
R •<— array of m empty lists 

i?[0] ^ MgnTrivalentGraphs{g ,n) first item contains all 3—valent graphs 
for fc in 1, \ldots, m — 1: 
Initialize R[k] as an empty list 
for G in R[k - 1]: 
for e in edge-orbits{G): 
if e is a loop: 

continue with next e 
G' ■(— contract{G, e) 
if G' not already in R[k]: 
append G' to R[k] 
return R 



Transforming boundary cycles under an isomorphism. The function 
transform_boundary_cycle is used when comparing marked fatgraphs: as 
the marking is a function on the boundary cycles, we need to know exactly 
which boundary cycle of the target graph corresponds to a given boundary 
cycle in the source graph. 

Recall that BoundaryCycle instances are defined as list of corners; func- 
tion transform_boundary_cycle takes a BoundaryCycle b and returns a new 
BoundaryCycle object b', obtained by transforming each corner according to 
a graph isomorphism. Indeed, transform_boundary_cycle is straightforward 
loop over the corners making up b: For each corner a new one is con- 

structed by transforming the vertex according to map pv, and displacing indices 
i and j by the rotation amount indicated by rot[\-f v] (modulo the number of 
edges attached to v). 

5 Generation of fatgraphs 

Let MgnGraphs be the function which, given two integers g, n as input, returns 
the collection of graphs. Let us fm'ther stipulate that the output result will 
be represented as a list R: the 0-th item in this list is the list of graphs with 
the maximal number m of edges; the fc-th item R[k] is the list of graphs having 
m — k edges. There are algorithmic advantages in this subdivision, which are 
explained below. 

Graphs with the maximal number of edges are trivalent graphs; they are com- 
puted by a separate function MgnTrivalentGraphs, described in Section [01 

We can then proceed to generate all graphs in 5lg,n by contraction of regular 
edges: through contracting one edge in trivalent graphs we get the list R[l] 
of all graphs with m — 1 edges; contracting one edge of G G wc get 

G' £ R[2] with TO — 2 edges, and so on. Pseudo-code for MgnGraphs is shown 
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in Algorithm [HI The loop at lines [51-[T^ is the core of the function: contract 
edges of the fat graph G (with m — k + 1 edges) to generate new fat graphs 
with m — k edges. However, we need not contract every edge of a fatgraph: if 
a E Aut G is an automorphism and e G E{G) is an edge, then the contracted 
graphs G" = G/e and G" = G/a{e) are isomorphic. Hence, we can restrict the 
computation to only one representative edge per orbit of the action induced 
by AutG on the set E{G); the edge.orbits function referenced at line |H] should 
return a list of representative edges, one per each orbit of Aut(G) on E{G). 

Lines [T^13l add G' to R[k] only if it is not already there. This is the most 
computationally expensive part of the MgnGraphs function: wc need to per- 
form a comparison between G' and each element in R[k]; testing equality of two 
fatgraphs requires computing if there are isomorphisms between the two, which 
can only be done by attempting enumeration of such isomorphisms. Fatgraph 
isomorphism is discussed in detail in Section 2] 

If Nk is the number of elements in R[k] and Tiso is the average time needed to 
determine if two graphs are isomorphic, then evaluating whether G' is already 
contained in R[k] takes 0{Nk -Tiso) time: thus, the subdivision of the output R 
into lists, each one holding graphs with a specific number of edges, reduces the 
number of fatgraph comparisons done in the innermost loop of MgnGraphs, 
resulting in a substantial shortening of the total running time. 

Note that the top-level function MgnGraphs is quite independent of the 
actual implementation of the Fatgraph type of objects: all is needed here, is 
that we have methods for enumerating edges of a Fatgraph object, contracting 
an edge, and testing two graphs for isomorphism. 

Lemma 5.1. If MgnTrivalentGraphs{g, n) returns the complete list of triva- 
Icnt fatgraphs in Jlg^n, then the function MgnGraphs defined above returns the 
complete set of fatgraphs ^g,n- 

Proof. By the above dissection of the algorithm, all we need to prove is that 
any fatgraph in 3ig,n can be obtained by a chain of edge contractions from a 
trivalent fatgraph. This follows immediately from the fact that any fatgraph 
vertex v of valence z ^ 3 can be expanded (in several ways) into vertices ui, V2 
of valences zi, Z2 such that z = {zi — 1) + {z2 — 1), plus a connecting edge. □ 

5.1 Generation of Trivalent Fatgraphs 

Generation of trivalent graphs can be tackled by an inductive procedure: given 
a trivalent graph, a new edge is added, which joins the midpoints of two existing 
edges. In order to determine which graphs should be input to this "edge addi- 
tion" procedure, one can follow the reverse route, and ascertain how a trivalent 
graph is transformed by deletion of an edge. 

Throughout this section, / and m stand for the number of vertices and edges 
of a graph; it will be clear from the context, which exact graph they are invari- 
ants of. 

5.1.1 Removal of edges 

Let G G 3^3, n be a connected trivalent graph. Each edge x G E{G) falls into one 
of the following categories: 
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Figure 8: Graphical illustration of fatgraph edge removal. Top row: a regular 
edge (crossed) is removed from an D^o,4 graph; its endpoints are further removed; 
the remaining edges are joined and the resulting graph is a trivalent fatgraph 
in !Ko.3. Bottom row: a loop is removed from a trivalent 3^o,4 graph; the stem 
together with its endpoints has to be removed as well; the remaining edges are 
joined, and we end up with a trivalent fatgraph in 3?o,3- 

A) X is a loop: both endpoints of x are attached to a single vertex v, another 
edge x' joins v with a distinct vertex w'; 

B) X joins two distinct vertices v,v' G V{G) and separates two distinct bound- 
ary cycles 13,13' G B(G); 

C) X joins two distinct vertices v,v' G V{G) but belongs to only one boundary 
cycle (3 G B{G), within which it occurs twice (once for each orientation). 

Deletion of edge x requires different adjustments in order to get a trivalent graph 
again in each of the three cases above; it also yields a different result in each 
case. 

Case A): If a; is a loop attached to v, then, after deletion of x, one needs to 
also delete the loose edge x' and the vertex v' (that is, join the two other edges 
attached to v'; see Figure |S1 bottom row). The resulting fatgraph G' has: 

> two vertices less than G: v and v' have been deleted; 

> three edges less: x, x' have been deleted and two other edges merged into 
one; 

> one boundary cycle less: the boundary cycle totally bounded by x has 
been removed. 

Therefore: 

2 - 2g' = x(G') = 1' -m' + n' 

= {l-2)- (m - 3) + (n - 1) 
^l^m + n^ X{G) =2-2(7, 
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hence g = g' , and 



(A) 



In case B), x joins distinct vertices v, v' and separates distinct boundary 
cycles (see Figure [51 top row) . Delete x and merge the two edges attached to 
each of the two vertices v and v'; in the process, the two boundary cycles /3, /3' 
also merge into one. The resulting fatgraph G' is connected. Indeed, given any 
two vertices u, u' S V{G'), there is a path (xi, . . . ,Xk) connecting u with u' in 
G. If this path passes through x, one can replace the occurrence of x with the 
perimeter — excluding x — of one of the two boundary cycles /3, /3' to get a path 
joining v and v' which avoids x, and thus projects to a path in G". Again we 
see that G" has: 

» two vertices less than G: v and v' have been deleted; 

» three edges less: x has been deleted and four other edges merged into two, 
pair by pair; 

» one boundary cycle less: the boundary cycles (3, (3' have been merged into 
one. 

Therefore g ~ g' , and 

G'e3l<,,„-i. (B) 

In case Cj, x joins distinct vertices w, v' but belongs into one boundary cycle 
(3 G B{G) only. Delete edge x and the two vertices v, v', joining the attached 
edges two by two as in case B). We distinguish two cases, depending on whether 
the resulting fatgraph is connected. 

C) If the resulting fatgraph G' is connected, then /3 e B{G) has been split into 
two distinct boundary cycles I3',f3" E B{G'). Indeed, write the boundary 
cycle /3 as an ordered sequence of oriented edges: yi ... — > 

Uk ^ yo- Assume the appear in this sequence in the exact order they 
are encountered when walking along /? in the sense given by the fatgraph 
orientation. The oriented edges are pairwise distinct: if yi and yj share 
the same supporting edge, then yi and yj have opposite orientations. By 
the initial assumption of case C), edge x must appear twice in the list: if x 
and X denote the two orientations of x, then yi = x and yj ~ x. Deleting 
X from j3 is (from a homotopy point of view) the same as replacing yi = x 
with X X, and yj = x with x x when walking a boundary cycle. 
Then we see that /3 splits into two disjoint cycles: 

/?' = yo ^ 2/1 ^ • ■ • ^ iji-i X ^ x^ yj+i >yk^yo, 

P" = yi+i > yj-1 ^ X -> S -> j/j+i. 

In this case, G' has: 

» two vertices less than G: v and v' have been deleted; 

» three edges less: x has been deleted and four other edges merged into 
two, pair by pair; 

» one boundary cycle more: the boundary cycle /? has been split in the 
pair p". 
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Therefore g' = g — 1 and n' = n + 1, so: 

G' G 3lg-i^n+l- (C') 

C") G' is a disconnected union of fatgraphs G'^ and Gj; for this statement 
to hold unconditionally, we temporarily allow a single circle into the set 
of connected fatgraphs (consider it a fatgraph with one closed edge and 
no vertices) as the one and only element of lRo.2- As will be shown in 
Lemma [5.21 this is irrelevant for the MgnTrivalentGraphs algorithm. 
Now: 

l[ + I2 = I — 2, m'^ + m'2 = m — 3, n[ + = n + I, 

hence: 

(2 - 2g[) + (2 - 2.9^) = - 2) - (m - 3) + (71 + 1) 
= {l-m + n) + 2 = 4:-2g 

So that g[ + 93 = g + 2, n[ + n'2 = n + 1, andjfl 

G' =G[®G'2e3ig'^,n[®^g',,n',- (C") 

5.1.2 Inverse construction 

If a; e E{G) is an edge of a fatgraph G, denote x and x the two opposite 
orientations of x. 

In the following, let 3?^ „ be the set of fatgraphs with a selected oriented 
edge: 

-.^ {{G,x) : G e Jlg^n,x e L{G)}. 
Similarly, let "Jig „ be the set of fatgraphs with two chosen oriented edges: 

3i;',„ :={(G,i,y) :G€3?g,„,x,y€L(G)}. 

The following abbreviations are convenient: 

Jl = UJlg^n, 3^' = „, Jl" = UJlg j^. 

Define the attachment of a new edge to a fatgraph in the following way. 
Given a fatgraph G and an oriented edge a;, we can create a new trivalent 
vertex v in the midpoint of x, and attach a new edge to it, in such a way that 
the two halves of x appear, in the cyclic order at v, in the same order induced 
by the orientation of x. Figure [3] depicts the process. 

We can now define maps that invert the constructions A), B), C) and C") 
defined in the previous section. 

Let pg n : — ^ ^g,n be the map that creates a fatgraph p{G, x) from a 
pair (G, x) by attaching the loose end of a "slip knot"F"l to the midpoint of x. 
The map p : D?' — >■ 3^ defined by p\ji> ^ := pg.„ is ostensibly inverse to A). 

To invert B) and C), define a map q : H" "R that operates as follows: 

®Here we use ® to indicate juxtaposition of graphs: Gi ® G2 is the (non-connected) fatgraph 
having two connected components Gi and G2. 

^"^A single 3-valent vertex with one loop attached and a regular edge with one loose end. 
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Figure 9: When adding a new vertex in the middle of an edge the cychc order 
depends on the oriented edge: the two orientations x and x get two inequivalent 
cychc orders. 

» Given (G, x, y) with x ^ y, the map q attaches a new edge to the midpoints 
of X and y; again the cychc order on the new midpoint vertices is chosen 
such that the two halves of x and y appear in the order induced by the 
orientations x, y. 

» When X = y, let us further stipulate that the construction of q{G,x,x) 
happens in two steps: 

(1) SL new trivalent vertex is created in the midpoint of x G E{G) and a 
new edge ^ is attached to it, 

(2) create a new trivalent vertex in the middle of the half-edge which 
comes first in the ordering induced by the orientation x; attach the 
loose end of the new edge ^ to this new vertex. 

It is clear that the above steps give an unambiguous definition of q in 
all cases where x and y are orientations of the same edge of G, that is, 
{G,x,x), {G,x,x), {G,x,x), and (G,x,x). 

Ostensibly, q inverts the edge removal in cases B) and C): the former applies 
when a graph G G 3?g,„ is sent to q{G) G Jig^n+i, the latter when G G IRg,„ is 
sent to q{G) G ^g+i^n-i- 

Finally, to invert C"), let us define 

From {G' ,x' ,G" ,x"), construct a new fatgraph by bridging G' and G" with a 
new edge, whose endpoints are in the midpoints of x' and x"; again, stipulate 
that the cyclic order on the new vertices is chosen such that the two halves of 
x' , x" appear in the order induced by the orientations x', x" . 
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Figure 10: Graphical illustration of maps p, q, Vg^n- Top left: p{G,x) attaches a "slipknot" to edge x. Top right: r2,5(G'i, x, 6*2, y) joins 
fatgraphs Gi and G2 with a new edge. Bottom: q(G,x,y) (left) and piG,x,y) (right); it is shown how changing the orientation of an 
edge can lead to different results. 



Summing up, any fatgraph G £ '^g,n belongs to the image of one of the above 
maps p, q, and r. There is considerable overlap among the different image sets: 
in fact, one can prove that r is superfluous. 

Lemma 5.2. Any fatgraph obtained by inverting construction C") lies in the 
image of maps p and q. 

Proof. Assume, on the contrarily, that G lies in the image of r only. Then, 
deletion of any edge x from G yields a disconnected graph G" ® G" . Both 
subgraphs G' and G" enjoy the same property, namely, that deletion of any 
edge disconnects: otherwise, if the removal of y £ E{G') docs not disconnect 
G", then neither docs it disconnect G = rg^n{G' ,G"), contrary to the initial 
assumption. As long as G' or G" has more than 3 edges, we can delete another 
edge; by recursively repeating the process, we end up with a fatgraph G* with 
Z* ^ 3 edges, which is again disconnected by removal of any edge. Since G* is 
trivalent, 3 • m* = 2-1*, therefore G* must have exactly 3 edges and 2 vertices. 
But all such fatgraphs belong to 3lo,3 or D?i_i, and it is readily checked that 
there is no way to add an edge such that the required property holds, that any 
deletion disconnects. □ 

5.1.3 The MgnTrivalentGraphs algorithm 

The stage is now set for implementing the recursive generation of trivalent 
graphs. Pseudo-code is listed in Algorithm [T] 

Lemma 5.3. MgnTrivalentGraphs{g,n) generates all trivalent fatgraphs for 
each given g, n. Only one representative per isomorphism class is returned. 

Proof. The function call MgnTrivalentGraphs{g,n) recursively calls itself to 
enumerate trivalent graphs of 3?(,_„„i and CHg^i^n+i- In particular, MgnTriva- 
lentGraphs must: 

» provide the full set of fatgraphs 3^0,3 and as induction base. 

» return the empty set when called with an invalid (g, n) pair; 

The general case is then quite straightforward: (1 ) apply maps p, q to every 
fatgraph in ^g^n-i, and q to every fatgraph in lRg_i^„_|_i; (2) discard all graphs 
that do not belong to 3?g,n; and (3) take only one graph per isomorphism class 
into the result set. 

To invert construction A), map p is applied to all fatgraphs G G 5lg^n-i', if 
a S Aut G, then p[a{G), a{x)) = p{G, x), therefore we can limit ourselves to one 
pair (G, x) per orbit of the automorphism group, saving a few computational 
cycles. Similarly, since g is a function of (G,x,y), which is by construction 
invariant under Aut G, we can again restrict to considering only one (G, x, y) 
per AutG-orbit; this is computed by the cdgc_pair_orbits{G) function. □ 

Note that there is no way to tell from G if fatgraphs p{G, x) and q{G, x, y) 
belong to ^g,n' one needs to check g and n before adding the resulting fatgraph 
to the result set R. 

The selection of only one representative fatgraph per isomorphism class can 
be done by removing duplicates from the collection of generated graphs in the 
end, or by running the isomorphism test before adding each graph to the working 
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Algorithm 7 Return a list of all connected trivalent fatgraphs with prescribed 
genus g and number of boundary cycles n. A fatgraph if "admissible" iff it has 
the prescribed genus g and number of boundary cycles n. 

dei MgnTrivalentGraphs{g, n): 

avoid infinite recursion in later statements 
if n= or {g,n) < (0,3): 
return empty list 

Induction base: Mo,3 and Mi^i 
if (5, n) = (0,3): 

return list of fatgraphs in 3lo,3 
ehf (5,n) = (1,1): 

return list of fatgraphs in Jli^i 

general case 
else: 
R <— empty list 

case A): hang a circle to all edges of graphs in 
for G in MgnTrivalentGraphs{g , n — 1): 
for X in edge-orbits{G): 
add p{G,x) to R if admissible 
add p{G,x) to R if admissible 

case B): bridge all edges of a single graph m Mg^n-l 
for G in MgnTrivalentGraphs{g , n — 1): 
for {x,y) in edge-pair.orbits{G): 
add q(G,x, y) to R if admissible 
add q{G,x, y)) to R if admissible 
add q{G, X, y) to R if admissible 
add (?(G, X, y) to R if admissible 

case C): bridge all edges of a single graph in Afg_i.„_|_i 
for G in MgnTnvalentGraphs{g — 1, n + 1): 
for {x,y) in edge-pair.orbits{G): 
add q(G,x, y) to R if admissible 
add q{G,x,y) to R if admissible 
add g(G, x, y) to R if admissible 
add q{G^ x, y) to R if admissible 

remove isomorphs from R 
return _R 
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list R. The computational complexity is quadratic in the number of generated 
graphs in both cases, but the latter option requires less memory. In any case, 
this isomorphism test is the most computationally intensive part of MgnTriva- 
lentGraphs. 

For an expanded discussion of the size of the result set R, and a comparison 
with other generation algorithms, see Appendix [Xj It would be interesting to 
re-implement the trivalent generation algorithm using the technique outlined in 
[23| . and compare it with the current (rather naive) algorithm. 

5.1.4 Implementing maps p{G,x) and q{G,x,y) 

Implementation of both functions is straightforward and pseudo-code is there- 
fore omitted^ the only question is how to represent the "oriented edges" that 
appear in the signature of maps p and q. 

In both p and q, the oriented edge x or ^ is used to determine how to 
attach a new edge to the midpoint of the target (unoricnted) edge x. We can 
thus represent an oriented edge as a pair (e, s) formed by a Fatgraph edge e 
and a "side" s: valid values for s are -f 1 and —1, interpreted as follows. The 
parameter s controls which of the two inequivalent cyclic orders the new trivalent 
vertex will be given. Let a, 6, c be the edges attached to the new vertex in 
the middle of e, where a,b are the two halves of e. If s is -1-1, then the new 
trivalent vertex will have the cyclic order a < 6 < c < a; if si is —1, then the 
edges a and b are swapped and the new trivalent vertex gets the cyclic order 
b < a < c < b instead. 

6 The homology complex of marked fatgraphs 

Bctti numbers of a complex (W*, can be reckoned (via a little linear algebra) 
from the matrix form D^^''> of the boundary operators Dk- Indeed, given that 
bk := dimHkiW,D) andHk{W,D) Zk{W, D)/ Bk{W, D) = Kci D^/ Dk-i{Wk-i), 
by the rank-nullity theorem we have dimKerl?/,; = dimVFfc — rankD^'^^ hence 
bk = dimKerDfc - dim Dk-i{Wk-i) = dimWk - ra-nk D^''^ - rankD^'^-^). 

In order to compute the matrix D^''\ we need to compute the coordinate 
vector of DkX^p for all vectors x^p in a basis of Wk- If (M^*, D*) is the fatgraph 

(k) 

complex, then the basis vectors x)j arc marked fatgraphs with k edges, and the 
differential Dk is defined as an alternating sum of edge contractions. Therefore, 
in order to compute the coordinate vector of DkX^ , one has to find the unique 

fatgraph x^'y^ which is isomorphic to a given contraction of x^*^^ and score a 
±1 coefficient depending on whether orientations agree or not. 

Although this approach works perfectly, it is practically inefficient. Indeed, 
lookups into the basis set {x'"^'^^ ^} of Wk-i require on average 0{N'^) iso- 
morphism checks. Still, we can take a shortcut; if two topological fatgraphs G 
and G' are not isomorphic, so are any two marked fatgraphs (G, i^) and (G", v'). 
Indeed, rearrange the rows and columns of the boundary operator matrix iD^^) 
so that marked fatgraphs (G, ly) over the same topological fatgraph G corre- 
spond to a block of consecutive indices. Then there is a rectangular portion 

^^The interested reader is referred to the publicly-available eodc at 
|http: //f atghol .googlecode ■ com] for details. 
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of D'^'^) that is uniquely determined by a pair of topological fatgraphs G and 
G' . The main function for computing the boundary operator matrix can thus 
loop over pairs of topological fatgraphs, and delegate computing the each rect- 
angular block to specialized code. There are n\/ \ Aut G\ marked fatgraphs per 
given topological fatgraph G, so this approach can cut running time down by 
0((n!)2). 

The generation of inequivalent marked fatgraphs (over the same topological 
fatgraph G) can be reduced to the (computationally easier) combinatorial prob- 
lem of finding cosets of a subgroup of the symmetric group ©„. In addition, the 
list of isomorphisms between G and G' can be cached and re-used for comparing 
all pairs of marked fatgraphs (G, v)^ (G', v'). This strategy is implemented by 
two linked algorithms: 

(1 ) MarkedFatgraphPool: Generate all inequivalent markings of a given 
topological fatgraph G. 

(2) compute_block: Given topological fatgraphs G and G', compute the 
rectangular block of a boundary operator matrix whose entries correspond 
to coordinates of D{G,v) w.r.t. {G\v'). 

6.1 Generation of inequivalent marked fatgraphs 

For any marked fatgraph (G, z^), denote [G, v] its isomorphism class; recall that 
B{—) is the functor associating a fatgraph with the set of its boundary cycles. 
Let N{G) be the sets of all markings over G and N{G) the set of isomorphism 
classes thereof: 

N{G) :={(G,z.) ^: B(G)^{l,...,n}}, 
N{G) :={[G,H S(G)^{l,...,n}}. 

Let (G, i>) be a chosen marked fatgraph. Define a group homomorphism: 

$:Aut(G)3ai — > D o B{a) o p-^ e (6) 

The set P ^ (f>(Aut G) is a subgroup of 6„. 

Lemma 6.1. The marked fatgraphs {G,P) and {G,aP) are isomorphic if and 
only if a E P. 

Proof. Let a <E P, then a^^ E P and there exists a G Aut G such that: 

a^^ = V o B{a) o , 

whence: 

(a o £/) o B{a) o ~ id, 

therefore a induces a marked fatgraph isomorphism between (G, D) and (G, crop). 

Conversely, let i> = and assume (G, P) and (G, v) arc isomorphic as 
marked fatgraphs: then there exists a G AutG such that P o B{a) o p"^ is the 
identity. Given any P o B{a') o G P we have: 

P 3 P o B{a') o p-i = P o (P^i o P) o B{a) o B(a)"^ o B{a') o P^^ 
= (p o p~^) o p o i?(a) o (P^^ o p) o B{a~^) o B{a') o p~^ 
= (7^^ o (p o -B(a) o P~^) o p o B{a^^ o a') o p~^ 

= cr"^ o (P o B{a^^ o a') o p-i) G ct^^P, 
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therefore P = a ^P, so a E P. 



□ 



Define a transitive action of 6„ over N{G) by cr • (G, i^) := (G, <t;^); this 
descends to a transitive action of S„ on N{G). By the previous Lemma, P is 
the stabiHzer of [G, i^] in N{G). 

Lemma 6.2. The action of 6„ on N{G) induces a bijective correspondence 
between isomorphism classes of marked fatgraphs and cosets of P in S„. 

Proof. Given isomorphic marked fatgraphs {G,v) and (G, i^'), let cr, cr' G (5„ 
be such that v = a o D and v' = g' o v. By definition of marked fatgraph 
isomorphism, there is a G Aut G such that the following diagram commutes: 



B{G) 



B{a) 



B(G) 



Hence commutativity of another diagram follows 

B{a) 



B{G) 



B{G) 



Thus (G, P) is isomorphic to (G, a~^a'oi?); therefore a^^a' G P, hence, cr' G crP, 
i.e., cr and cr' belong into the same coset of P. 
Conversely, let r, r' G aP; explicitly: 

T = a o D o B{a) o D^^ , t' — a o D o B{a') o . 

Set u = T o D, v' ~ t' o D] substituting back the definition of r, we have: 

v ^ a o D o B{a) o D^^ o v = a o u o B{a), 

whence D = cr^^ o j/ o _B(a)~^, and: 

ly' ^ (7 o p o B{a') = a o (cr^^ o i/ o _B(a)^"'^) o B{a') = v o B{a^^ o a'), 

therefore o a' is an isomorphism between the marked fatgraphs (G,v) and 
{G,u'). □ 

The following is an easy corollary of the transitivity of the action of 6„ on 
N{G). 

Lemma 6.3. Given any marking v on the fatgraph G, there exist a G ©„ and 
a G Aut G such that: v = a o V o a. 

Proof. By Lemma [6.21 there exists a G such that [G^v] = [G, ct o V], i.e., 
(G, v) is isomorphic to (G, cr o z/) . If a G Aut G is this fatgraph isomorphism. 
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then the following diagram commutes: 



B{G) ^ B{G) 



Therefore i> = a o D o B{a). □ 

The MarkedFatgraphPool algorithm. Given a fatgraph G and a Fatgraph 
object G representing it, let us stipulate that D be the marking on G that enu- 
merates boundary cycles of G in the order they are returned by the function 
compute-boundary.cycles{G) . By Lemma lOl every (G, z^^-'^) can then be ex- 
pressed (up to isomorphism) as (G,cr^^-' op) with a^^^ G S„. The set {cr^^'} enu- 
merates all distinct isomorphism classes of marked fatgraphs over G iff {cr'-'^P} 
runs over all distinct coscts of P in 6„ (by Lemma 16. 2p . 

The MarkedFatgraphPool function computes the set N{G) of isomorphism 
classes [G, ly] . 

Theorem 6.1. Given a Fatgraph G as input, the output of MarkedFatgraphPool{G) , 
as computed by Algorithm \^ is a tuple {graph, P, A, markings, oricntable) , 
whose components are defined as follows: 

» The graph item is the underlying Fatgraph object G. 

» The P slot holds a list of all elements in the group P ~ $(AutG). 

» A corresponding set of pre-image representatives ( each element is an au- 
tomorphism of G) is stored into A: permutation P[i] is induced by auto- 
morphism A[i], i.e., if TT ~ P[i] and a = A[i] then n = ^{a). 

> The markings item holds the list {cr^^'} of distinct cosets of P (representing 
inequivalent markings). 

» orientable is a boolean value indicating whether any (G, v) in the pool is 
orientableV^ 

We need a separate boolean variable to record the orientability of the family 
of marked fatgraphs N{G) = {[G,v)}, because the automorphism group of a 
marked fatgraph Aut(G, i^) can be a proper subgroup of AutG: hence, {G,v) 
can be orientable even if G is not. 

Proof. Generation of all inequivalent markings over G is a direct application of 
Lemma 16.21 performed in two steps: 

(1 ) In the first step, for each automorphism a G Aut G compute the permu- 
tation $(a) it induces on the set of boundary components, and form the 
subgroup P. The subgroup P and the associated set of automorphisms 
A C Aut G are stored in variables P and A. 

^^It is an immediate corollary of Lemma 16.31 that if one marked fatgraph (G, t'*) has an 
orientation-reversing automorphisms, then every marked fatgraph {G, u) over the same topo- 
logical fatgraph G has an orientation-reversing automorphism. 



41 



Algorithm 8 Compute the distinct markings of a given fatgraph. Input to the 
algorithm is a Fatgraph object G; final result is a tuple (G, P, A, markings, 
orientable) which represents the set N{G) of isomorphism classes of marked 
fat graphs. 

def phi{a, G): 

TT <— array of n elements 

for srcdndex, src^cycle in enumcratc{G.boundary_cycles): 
dst-cyclc <— a.transform.boundary.cycle{src-cycle) 
if dst_cycle not in G .boundary .cycles: 

abort and signal error to caller 
else: 

dstJndex <— index of dst-cycle in G.boundary_cycles 
TT[srcJndex] ^ dstJndex 
return tt 

def MarkedFatgraphPool{G): 

F <— empty list 
A ■(— empty list 

assume (G, v) is orientable until we have counter— evidence 
orientable True 

step (1): loop over AutG 
for a in G.automorphisms{): 
try: 

TT ^ plii{a, G) 
except phi failed: 

continue with next a 
if permutation tt is identity: 
found a new automorphism: 

— does it reverse orientation? 
if a.is-oricntation^cversing{): 

orientable ^ False 

— does it define a new marking? 
if TT not in P: 

append tt to P 
append a to A 
step (2): enumerate cosets of P 
markings <r- [ ] 
for a in ©„: 
for TT in P: 
if TT o (7 in markings: 

continue with next a 
add a to markings 
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(2) In the second step, compute cosets of P by exhaustive enumeration. They 
are recast into the list {cr'^^^}, which is stored into the markings variable. 

As an important by-product of the computation, the automorphism group Aut(G', D) 
is computed, and used to determine if the marked fatgraphs in the pool arc ori- 
cntablc. 

The auxiliary function phi computes the permutation $(a) = Do B{a) oi>~^. 
A permutation tt is created and returned; it is represented by an array with n 
slots, which is initially empty and is then stepwise constructed by iterating over 
boundary cycles. Indeed, the boundary cycle src-cycle is transformed accord- 
ing to B{a) and its position in the list of boundary cycles of G is then looked 
up. Note that this lookup may fail: there are in fact cases, in which the Fat- 
graph, isomorphisms algorithm finds a valid mapping, that however does not 
preserve the markings on boundary cycles; such failures need to be dealt with 
by rejecting a as a Fatgraph automorphism. 

Step (1 ) of the computation is performed in lines [T5H?7l 

» Computation of the permutation tt (induced by a on the boundary cycles 
of G) may fail; if this happens, the algorithm ignores a and proceeds with 
another automorphism. 

» If a preserves the boundary cycles pointwise, then it induces an automor- 
phism of the marked graph and we need to test whether it preserves or 
reverses orientation. 

» There arc |Ker$| distinct automorphisms inducing the same permutation 
on boundary cycles: if tt is already in P, discard it and continue with the 
next a. 

By Lemmain^l there are as many distinct markings as there are cosets of P in 
©„. Step (2) of the algorithm proceeds by simply enumerating all permutations 
in &n, with marking initially set to the empty list; for each permutation a a 
test is made as to whether aP intersects the list markings (lines [351 - I37)) : if it 
does not, then the marking induced by a is added to the list. 

□ 

A constructive version of Lemma l6.3l can now be implemented: the following 
function index_and_aut, given a Fatgrapli object G and a marking, returns the 
permutation (by index number j in G. markings) and fatgraph automorphism 
a = G.A[i] such that the topological fatgraph G decorated with marking is 
isomorphic (through a) to the same graph decorated with G.markings[j]. 

def index_and_aut{G, marking): 
for (i,7r) in enumerate{G.P): 

r (7 O TT 

if T in G. markings: 

j <— index of r in G. markings 

return (j, 
else: 

continue with next n 

The algorithm enumerates all permutations tt € P, and compares cr o tt to every 
clement of G. markings: by Lemma 16.21 we know that one must match. 
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6.2 Computing boundary operator matrix blocks 

The differential D{G, v) is computed by summing contractions of regular edges 
in G (with alternating signs); likewise, the matrix block corresponding to co- 
ordinates of the families of marked fatgraphs {(G, i^)} and {{G' ,v'Y\ can be 
decomposed into a sum of blocks, each block representing the coordinates of 
{(G/e, i^)}e6_B(G) projected on the linear span of {(G', v')}- 

More precisely, given any two fatgraphs Gi (with m edges) and G2 (with 
m — 1 edges), let Xi,X2 C ^g_„ be the linear span of N{Gi) and N{G2) 
respectively, and denote by prj^^^ the linear projection on subspace X2. Recall 
that, for any fatgraph G, we have D{G) = ^ ±d''^^(G), where the sum is taken 
over all regular edges e of G, and d'^^ is the contraction of edge e. 

Let G be the fatgraph obtained by contracting the chosen edge e in Gi. If 
G2 and G are isomorphic, then the three graphs are related by the following 
diagram of fatgraph morphisms, where /i is the contraction map and /2 is a 
fatgraph isomorphism: 

Gi (7) 



G — G2 

The above diagram ([7]) functorially induces a diagram on the set of boundary 
cycles: 

B(Gi) (8) 

B{fi 

B{G) ^B(G2) 

Diagram commutes iff /i, /2 can be extended to morphisms of marked 
fatgraphs /i : {Gi,vi) (G,v) and h ■ {G,v) {G2,V2)- 

Now choose Fatgraph objects Gi, G, G2 representing Gi, G, G2. 

Let i^i, P, V2 be the markings on Gi, G, G2 that enumerate boundary cycles 
in the order they are returned by the function computejDOundary_cycle applied 
to Gi, G, G2 respectively. Define 01,02 £ ©n by: 

01 := o -B(/i) o \ 02 := o -B(/2) o (9) 

Lemma 6.4. Given any marking vi on Gi, choose ai G 6„ such that vi =^ ui o Ui 
and define: 

z^2 := CTi o 0~^ o 0~^ o ^2, (10) 
Then V2 is the unique marking on G2 such that diagram ([8]) commutes. 

Proof. Let a2 '.= (Ji ° 4>i^ ° 4>2^ ■ We need to prove that the external square in 
diagram is commutative; indeed, we have: 

172 = CTi o {vi o B{fi)~^ o jy"i) o{Do B{f2y^ o i?2) 

= aivi o B{f2 o fiy^ o p2"\ 
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so that: 



V2 ° B{f2) O B{fi) = (72 O 7^2 O B[f2 O /i) 

= aivi o B{f2 o fiy^ o zy~^ o 1/2 o B{f2 o /i) 

= (Tl O £/l = i^l . 

The uniqueness assertion is of immediate proof, since maps -B(/i) and B{f2) 
are invertible. □ 

Let pi, P2 be the MarkedFatgraphPool output corresponding to Gi, G2, 
and let {vi^}j=i^...^Nn {i^2'''}k=i,....N2 be the enumeration of fatgraph markings 
corresponding to items in the hsts pi. markings and p2.inarkmgs respectively. 

Lemma 6.5. For any regular edge e of Gi, and any choice of j 6 {1, . . . , Ni}, 
there exist unique k € {1, . . . , A''2} and s £ { — 1, 0, +1} such that: 

pr^^(d(^)[Gi,i^(^')]) =s-[G2,4'^]. (11) 

Proof. If G2 and G ~ Gi/e are not isomorphic, then, for any marking vi, 
d{Gi,vi) has no component in the subspace X2 = {(62,2^2)}, so the assertion 
is true with s = 0. 

Otherwise, by Lemma 16.41 given vi ~ there is a unique 1^2 such that s 
can be non-null; by Lemma 16.31 there exist v''^'' :~ ctj''^ o D2 and a G Aut G such 
that: 

(1) the marked fatgraph (G2, 1^2^^) is a representative of the isomorphism class 
[02,1^2]; 

(2) a gives the isomorphism between marked fatgraphs (G2, 1^2) and (G2, '^2'^''); 

(3) v^^ is the marking on G2 represented by fc-th item in list p2. markings. 

The coefficient s must then be ±1 since both (G2,i^f') and d^'^{Gx,v(^) are 
(isomorphic to) elements in the basis of X2. □ 

Theorem 6.2. Given MarkedFatgraphPool objects pi, p2, and a chosen edge 
e of Gi, the function compute-block in Algorithm \^ returns the set S of all 
triplets {j, k, s) with s = ±1 such that: 

Wx.{d^''nGu>^['^])=s-[G2,4% (12) 

Proof. The algorithm closely follows the computation done before Lemmas 16.41 
and in the proof of Lemma 16.51 

If G2 and G = Gi/e are not isomorphic, then d'-'^)[Gi, i^i] has no component 
in the subspace X2 generated by {[G2,t^2]}, whatever the marking i/i. The 
assertion is thus satisfied by S* = 0, i.e., an empty list of triplets (j, k, s) (linesjS]- 
Elin Algorithm O. 

If G2 is isomorphic to G = Gi/e through /2, then Lemma [6.41 provides the 
explicit formula i^^'''' ~ a[^^ o (f)~^ o o D2, where a^^'^ = v^^"^ o v^^ . 
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Algorithm 9 Return the set S of triplets (j, k, s) such that eq. ([T^ holds for 
{Gi,h'[''^) and (02,1^2''^) obtained by contracting e in all marked graphs in pi 
and projecting onto graphs in the p2 family. 

def conipute_block{pi, e, P2): 
Gi ^ pi. graph 
G2 P2-graph 
G -h- contract{Gi, e) 
if G and G2 arc not isomorphic: 

return empty list 
else: 

result empty list 

f2 ■v- first isomorphism computed by Fatgraph.isomorphisms{G , G2) 
(pi^ -f" compute_philJnv{G , Gi, e) 
(j)2^ ^ computc-phi2Jnv{G, G2, f2) 
for (j, cr) in cnumcratc{pi. markings): 

{k, a) indcx_and_aut{p2, cr o 4>i^ o (j)^^) 

p ^ Gi.orkmt[e] 

s 4— (—1)^ * compare-orientations{f2) * compare-orientations{a) 
append (j, k, s) to result 
return result 

def compute_phil_inv{G, Gi, e): 

T -i— empty array of n elements 

for i, 6 in enumerate{Gi. boundary ^cycles): 

b' contract_boundary_cyclc{Gi, b, e) 

i' index of b' in G .boundary _cycles 

T[i'] <— i 
return r 

def compute_phi2_inv{G, G2, f2)' 
t' empty array of n elements 
for i, b in enumerate(G2-boundaryjcycles): 

b' -i— transform_boundary_cycle{f2, b) 

i' <— index of b' in G .boundary _cycles 

t'[i'] f- I 
return r' 
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By assumption, i/i numbers the boundary cycles on Gi in the order they are 
returned by running function compute_boundary_cycles on Gi, so is the 
permutation corresponding to the j-th element in pi. markings. 

The map 4>i is easy to compute: again, given that both D and Di number the 
boundary cycles of G and Gi in the order they are returned by computc_boundary_cycles, 
the auxiliary function computc_philJnv incrementally builds the result by loop- 
ing over Gi.boundary^cycles, contracting the target edge, and mapping the cor- 
responding indices. 

Computation of the map (f>2 depends on the isomorphism /2; however, two 
different choices for /2 will not change the outcome of the algorithm: in the final 
loop at lines [r^fTOl only the sign of /2 is used, and the sign is constant across 
all isomorphisms having the same source and target fatgraphs (iff they are both 
orientable). Computation of (in the auxiliary function compute_phi2Jnv) 
is done in the same way as the computation of (p^^ , except we transform b to b' 
by means of transform_boundary_cycle{f2, — ), i.e., B{f2). 

Finally, for every marking cr^' in pi. markings (representing j^^^), we know 

by Lemma 16.51 that there is a unique index k and a G Aut G2 such that: a[^'^ o 

(p^^ o (j)^^ = CTj*^^ is the fc-th item in p2.markings (representing 1^2'^'^), and such 
that the following chain: 

Gi G ^ G2 ^ G2 

fl h " 

extends to a marked fatgraph morphism: 

(Gi, z/p")) — ^ (G, v) ^ (G2, V2) ^ (G2, ^ 

fl f2 " 

The sign s is then obtained by comparing the orientation uj2 of (G2 , 1^2'^^ ) 
with the push-forward orientation (a o /2 o fi)^uji, where uti is the orientation 

on (Gi, i^P^), and multiplying by the alternating sign from the homology differ- 
ential. There are four components that make up s: 

» the sign given by the contraction /i: this is -1-1 by definition since the 
"child" fatgraph G inherits the orientation from the "parent" fatgraph 
Gi; 

» the sign given by the isomorphism /2: this is obtained by comparing (/2)*w 
with LU2 , which is implemented for a generic isomorphism by the function 
compare-orientations; 

» the sign of the automorphism a of G2 which transforms the push-forward 
marking into the chosen representative in the same orbit: this again can 
be computed by comparing (a)*aj2 with lj2 and only depends on the action 
of a on edges of G2 ; 

» the alternating sign from the homology differential, which only depends 
on the position p of edge e within the order wi . 

The product of the three non-trivial components is returned as the sign s 
(line [IS]). □ 
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Algorithm 10 Compute the boundary operator matrix, block by block. 



dcf conipute_boundary_operator{m, graphs): 
Ni ■<— number of graphs with m edges 
N2 ■<— number of graphs with m — 1 edges 
£)("') ^ Ni X N2 matrix, initially null 
jo ^ 

for Gi in graphs[m]: 
Pi ^ MarkedFatgraphPool{Gi) 

fco ^ 

for G2 in graphs[m~l]: 
P2 ^ MarkedFatgraphPool{G2) 
for e in Gi. edges: 
if e is a loop: 

continue with next e 
for {j,k,s) in coinpute-block{pi, e, P2): 
add s to entry D'™) [fc + fco , j + jo] 
increment kg by the number of inequivalent markings in p2 
increment jo by the number of inequivalent markings in pi 
return D^") 



6.3 Matrix form of the differential D 

The compute-boundary^operators function (Algorithm llOp computes the matrix 
form of the differential D restricted to the linear space generated by 

fat graphs with m edges. 

Input to the function are the number m and the list of graphs, divided by 
number of fatgraph edges: graphs[m\ is the list of fatgraphs with m edges. 

The output matrix Z)*^'") is constructed incrementally: it starts with all 
entries set to 0, and is then populated blockwise. Indeed, for every pair of 
MarkedFatgraphPool objects pi (from a graph with m edges) and p2 (with m — 1 
edges), and every non-loop edge e, the rectangular matrix block whose upper- left 
corner is at indices jo, fco is summed the block resulting from computc_block{pl, e, p2). 



7 Conclusions 

A Python implementatiorF^ of the algorithms presented in this paper has been 
actually used to compute the Betti numbers of all Mg^„ with 2g + n ^ 6. The 
results are summarized in Table [TJ Corresponding calculations based on the- 
oretical results are scattered across a wide array of publications. For g ^ 1, 
the groups H^{l\{g^n, Q) are known from the works of Mumford [^5] and Harer 
|17) : i?^(Mg_„,Q) has been computed also by Harer in [T7]; a comprehensive 
statement with a new proof is given by Arbarello and Cornalba in [5] (where 
a minor mistake in Harer's statement is corrected). The complete homology of 
Mi^2 and M2.1 has been published in Godin's paper [M]. The homology of the 

'^^Code publicly available at |http: //code .google ■ com/p/f atghol"| 
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No. of edges: 12 11 10 9 8 1 6 5 4 3 2 
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Table 2: Number of distinct abstract fatgraphs with the given genus g and 
number of boundary cycles n. For readability, null values have been omitted 
and the corresponding entry left blank. 



Mo,* spaces is computed in [TTl Corollary, 3.10]; see [T] for an alternative ap- 
proach using results from [28| to compute the Poincare polynomial of Mo,ti- The 
Poincare-Serre polynomial of M2.2 follows as a special case of Corollary III. 2. 2 
in Tommasi's [57]; the results also follows by combining [T^l p. 22] with [T51 
Appendix A]. The rational cohomology of Mi. 4 is completely described in The- 
orem 1 of [TS]; the Betti numbers were already present in [T31 p. 22]. In all these 
cases, the numerical results agree with the values in Table [TJ 

An internal verification step in the code computes the classical and virtual 
Euler characteristics of the fatgraph complex; the computed values match those 
published in [121 [HI [S], where they are derived by theoretical means. 

As a side effect of the computation, the entire family of fatgraphs Olg_n (with 
2g + n ^ 6) has been computed, and for each fatgraph the isomorphism group is 
known. The full list of fatgraphs and their isomorphisms is too long to print here. 



but the data is publicly available at http : //f atghol . googlecode . com/download/list 



Tables [5] and 13] provide a numerical summary of the results. 



7.1 Performance 

Table U] gives a summary of the running times obtained on the idhydra.uzh. ch 
cluster at the University of Zurich; Figure [TT] provides a graphical representation 
of the same data. The computational demands of the code are such that the 
homology of Mg.„ can actually be computed on desktop-class hardware for 2g + 
n < 6. 

The scatter plot in Figure [TT] shows that the time spent in computation 
of the ID^™) matrix ranks done in Stage III can become the dominant contri- 
bution to the total running time as the number of fatgraphs increases. This 
highlights a limitation of the program: the large number of fatgraphs in the 
Kontsevich complex might turn out to be a challenge for today's sparse linear 
algebra software. 

However, the set of fatgraphs for a given {g,n) pair has to be generated 
prior to computing the matrices D^™''^: a very large set of graphs can exhaust 
the computer's memory long before computation time becomes a blocking issue. 
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Total 
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327 


31262 


4583322 


2 


37 


4168 


747664 


142 


43054 



Table 3: Number of distinct orientable marked fatgraphs in the Penner- 
Kontsevich complex of each of the indicated Mg^„ spaces. For readability, null 
values have been omitted and the corresponding entry left blank. 



Time (s) 


: Stage I 


Stage II 


Stage III 


Total 


Mo,3 


< 1ms 


< 1ms 


0.03 


0.12 


Mo,4 


0.05 


0.09 


< 1ms 


0.29 


M0.5 


4.78 


21.91 


1.85 


29.43 




2542.56 


16011.70 


179157.39 


233007.06 


Mi.i 


< 1ms 


< 1ms 


0.010 


0.128 


M1.2 


0.05 


0.08 


< 1ms 


0.27 


Mi,3 


40.56 


136.88 


< 1ms 


174.75 


Mi,4 


82486.51 


336633.75 


4872.69 


424615.85 


M2^1 


2.39 


4.76 


< 1ms 


7.39 


M2,2 


43402.18 


181091.11 


5.57 


224694.61 



Table 4: Total CPU time (seconds) used by the Betti numbers computation for 
the indicated Mg^„ spaces. The C++ library LinBox [22j[9] was used for the rank 
computations in Stage III. Running time was sampled on the idhydra . uzh . ch 
computer of the University of Zurich, equipped with 480GB of RAM and Intel 
Xeon CPUs model X7542 running at 2.67GHz; Python version 2.6.0 installed 
on the SUSE Linux Enterprise Server 11 64-bits operating system was used to 
execute the program. The system timer has a resolution of 1ms, but times are 
less accurate than that, because of the scheduling jitter in multitasking systems. 
The "Total" column does not just report the sum of the three stages, but also 
accounts for the time the program spent in I/O and memory management. 
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stage I 
in □ Stage II 
J - o Stage III 
A Total 



M2,l 



~r~r 



142 327 4168 31262 

Number of marked fatgraphs 



747664 4583322 



Figure 11: Scatter plot of the data in Tabled Both axes use log-scale. Note how 
Stage III (computation of the boundary operators rank) becomes the dominant 
task as the number of marked fatgraphs increases. 

7.2 Application to other fat graph complexes 

In [14] , V. Godin defined a "bordered fatgraph complex" , which computes the 
integral homology of the moduli spaces of Riemann surfaces with boundaries. 
Godin's fatgraphs extend the abstract fatgraph by requiring that a leaf (i.e., 
a univalent vertex), and only one, is present in each boundary cycle. The 
bordered fatgraph complex is then constructed exactly as the fatgraph complex 
presented here, with the proviso that an edge ending in a univalent vertex is 
never contracted: hence, the differential D is given by the sum of contraction 
of non-loop non-leaf edges. 

The algorithms of this paper can easily be adapted to compute the homology 
of Godin's bordered fatgraph complex: after generating the family of marked 
fatgraphs of a given {g,n) pair, we decorate each marked fatgraph with leaves; 
compute the matrix form of the differential D and then reduce it to Smith 
normal form to reckon the integral homology modules of the moduli space of 
bordered surfaces. 

There is no need for checking duplicates in the set of bordered fatgraphs 
so generated^ therefore the decoration step can be implemented efficiently. A 
shortcut can also be taken in computing the matrix D: since leaf edges arc 
never contracted, the differential on bordered fatgraphs can be deduced easily 
from the differential on marked fatgraphs. However, the number of bordered 

^*If two "bordered fatgraphs" were isomorphic, they would remain such if we remove the 
leaves and the edge supporting them, which would give us isomorphic marked fatgraphs 
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fat graphs is much larger than the number of marked fatgraphs0this means that 
the final linear algebra computations require even more computational resources 
than they do for Mg^„ computations. 

7.3 Future development directions 

There are a number of directions in which the current algorithms and code could 
be improved. 

As already noted, the generation algorithms produce quite a number of du- 
plicates, that have to be removed using a quadratic-complexity procedure. A 
variant of the "isomorph-free generation" algorithm of McKay [13] could replace 
the naive MgnTrivalentGraphs code; the question of which algorithm would be 
faster has probably to be sorted out empirically, the critical performance factor 
being the number of times the "isomorphism" test is invoked. 

Another approach would be to turn the generation procedure "upside down" : 
instead of starting with trivalcnt graphs and contracting edges, one could start 
with (g, n)-fatgraphs with one vertex and expand those until the whole set of fat- 
graphs is generated. This would have the advantage that the chromatic fatgraph 
polynomial of BoUobas and Riordan [7] is available as an invariant to speed up 
the isomorphism procedure. On the other hand, the number of fatgraphs gener- 
ated this way seems consistently larger than the number of fatgraphs generated 
with the procedure adopted here (see Section [A. 2[) . 

So far, the major obstacle to applying the algorithms of this paper to a wider 
range of moduli spaces has been the large number of fatgraphs involved: it affects 
both the total run time and memory consumption of the code. Most algorithms 
described here lend themselves naturally to parallelization, so it would be pos- 
sible to rewrite the program to exploit several processors and distributed mem- 
ory, which could solve both issues. However, the number of generated fatgraphs 
grows supcr-exponentially in the asymptotic limit [31 0], so any implementa- 
tion of the algorithms outlined here will soon hit the limit of any present-day 
computing device. The question remains open, whether more significant result 
could be obtained before hitting the limits of today's computers. 
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leaf may be regarded as a choice of an edge or a vertex along a boundary cycle; if there 
are pi vertices (counted with multiplicities) and qi edges along the j-th boundary cycle, then 
the number of ways we could possibly add leaves to a marked fatgraph G is rir2 ■ ■ ■ Vn, with 
= Pi + Qi so that; 

ri -I- r-2 -I- ■ ■ ■ -I- r„ = ^ pi + ^ qt = ^ 2„ -|- 2m = Am, 

i i uev(G) 

where m is the total number of edges and is the valence of vertex v. 



52 



I am specially indebted, and would like to express my gratitude, to Enrico 
Arbarello, Domcnico Fiorenza and Gilberto Bini, for their constant encourage- 
ment and manifold support during the long preparation of this paper. 

I also would like to thank Orsola Tommasi and Dan Petersen, who pointed 
out inaccuracies in an earlier versions of this text and suggested improvements. 

A Comparison of fatgraph generation methods 

This section compares three different approaches to generating trivalent fat- 
graphs: namely, we compare the Mgn Trivalent Graphs algorithm described in 
Section fS . 1 .31 with two alternatives!^ Tabic [5] presents a summary of results. 

None of the suggested algorithms is capable of directly producing an isomorph- 
free set of distinct fatgraphs; they all produce a larger set of fatgraphs that 
must be reduced by taking only one representative per isomorphism class of fat- 
graphs. Therefore, Table [5] also reports the actual number of distinct fatgraphs 
for a given n pair; not all counts are known: a cell is left empty when the cor- 
responding count has not yet been computed. From the results gathered so far, 
it is apparent that all algorithms overestimate the actual number of fatgraphs. 

In what follows, let N{g,n) := |IRg.„| be the number of distinct {g,n)- 
fat graphs; also define: 



It is trivial to check that mmax and m^i^ are the maximum and minimum 
number of edges that a (g, ri,)-fatgraph can have. 

A.l Generation by recursive edge addition 

The algorithm MgnTrivalentGraphs described in Section [5 . 1 . 31 produces a (g, n)- 
fatgraph by adding an edge to fatgraphs with lower {g,n); the procedure can 
then be applied recursively. 

Let Ni(g, n) be the number of (non distinct) fatgraphs returned by MgnTrivalcntGraphs(g ,n) . 
According to Section [STjSl this can be written as: 



where Ni^a, -^i.c are the mmrbers of fatgraphs constructed in cases A), 

B), G') of Algorithm El 

In case A), wc have 1 generated (g, n)-fatgraph per each pair formed by a 
(g, n— l)-fatgraph and one of its oriented edges, modulo the action of the auto- 
morphism group Aut(G). However, the number of orbits of this Aut(G)-action is 
difficult to estimate. Since the generic fatgraph only has one automorphism, an 

^^The author is aware of no other algorithm for generating the set of all fatgraphs. The 
comparison here is taken with the solutions used in earlier attempts of implementation of the 
iFatGHoLi software. 



m^^^{g,n) 



2g + n, 

6,g + 3n — 6 = 3^ — 6, 
2g + n-l^i-l. 



(13) 



Ni[g, n) = Ni^Aig, n) + Ni^b{9, n) + 7Vi,c(<?, n), 
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Table 5: Number of (non- unique) trivalent fatgraphs generated according to 
different algorithms. The TV column reports the actual number of distinct fat- 
graphs for the given g, n; empty cells mean the corresponding number has not 
been computed. The column lists upper bounds for the recursive generation 
algorithm MgnTrivalentGraphs (see Listing[7l); values marked with the "f" sym- 
bol are estimated using earlier values of because the corresponding values of 
N are not available. The N2 values bound from above the number of fatgraphs 
generated by grafting binary trees into clovers. Finally, N3 is the count of fat- 
graphs generated by enumerating pairs of permutations (as per combinatorial 
definition of fatgraph) . 

upper bound can instead be given by considering all pairs formed by a fatgraph 
and an oriented edge: 

Ni^Aig, n) < iVi^^(g, n) := 2 • mn^^^{g, n - I) ■ N{g, n - 1). 

In case B), the algorithm generates one (5, n)-fatgraph per each triplet 
formed by a {g,n — l)-fatgraph and two oriented edges, not necessarily dis- 
tinct (modulo the action of AutG); a similar remark about the upper bound 
applies: 

Ni^Bia.n) ^ N+g{g,n) (2 • m„,ax(5," ^ 1))' ' ^^(.9," - 1) 

In case C), the computation is exactly the same, except we apply the q 
construction to fatgraphs belonging in 3ig^i_n+i' 

iVi,c(5, ") ^ N+cig, n) := 4 • m„,ax(5 - l,n + 1)^ ■ N{g - l,n + I). 
Table [S] shows the upper bound given by 

N+{g, n) N+^{g, n) + N+g{g, n) + N+^ig, n). 

According to Table [51 the MgnTrivalentGraphs algorithm outperforms the 
alternative procedures when 2g + n < 7, and apparently generates a much 
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larger set of fatgraphs when 2g + n > 7. However, the values were obtained 
usmg N^{g,ri) instead of N{g,ri) in recursive computations when the actual 
value of N{g, n) is not known; therefore N^{g, n) might grossly overestimate the 
number of graphs considered by MgnTrivalentGraphs for 2g + n > 6. Further 
investigation is needed to ascertain whether this is due to the algorithm of 
Section rA.2l being asymptotically faster, or to the estimate for Ni{g,n) being 
grossly imprecise when no data about the real number of trivalent fatgraphs in 
the recursion step is known. However, the author conjectures that this estimate 
holds: 

Ni{g,n)^0{e)-N{g,n) (14) 

A. 2 Generation by insertion of binary trees 

A different approach is the following: 

» Generate all distinct (g^ n)-fatgraphs with 1 vertex; each such fatgraph 
has TOmin(ff, n) edges, hence the vertex has valence 2 • m^in{g, n). 

» Given any such fatgraph Gq, build a trivalent (5, ri)-fatgraph G by replac- 
ing the vertex with a full binary tree on 2 • r7imin((7, n) leaves. 

Call a fatgraph with only one vertex a clover. Let N2{g,n) be the number of 
distinct (17, n)-clovers; we can estimate it as follows. 

Lemma A.l. The number of isomorphic clovers is equal to the number of orbits 
of the adjoint action 0/ (12 . . . 2m) over the set of self-conjugate permutations 
Wi e 62m : (rf = id}. 

Proof. Let Gq = (L; croi fi? 0'2) be a (g, ri)-fatgraph given in combinatorial form, 
where L = {1, . . . , 2m} and tXi G 6m. If Go is a clover, then ctq is a permutation 
formed by just one cycle; without loss of generality we may assume ctq is the 
rotation (12... 2m). Let Gq = (L; CTq, tri, (T2) be another (5, n)-clover: by the 
same reasoning we have ctq = ^0 = (12 . . . m); if / : Gq — Gq is an isomorphism, 
then / commutes with cfq hence f = for some j |2m. Therefore, from ci o / = 
/ o (j'l we get cTj = (Tq o (Tx o (Tq . This proves the claim. □ 

Lemma A. 2. Let L be a finite set of I = p ■ q elements. The number of 
permutations of L which can be expressed as product of q disjoint p- cycles is: 

C{p,q)=f[liipi-j). (15) 

1=1 j=l 

Proof. Without loss of generality we can assume L ~ {1, . . . ,pq}; let t G Gpq 
be a permutation composed of q disjoint p-cycles. We can give a "canonical" 
form to T if we order its cycles by stipulating that: 

» a cycle (aia2 . . . Op) is always written such that ai = minoi; 

» {aia2 . . . Up) precedes (6162 . ■ - bp) iff min cii < min bi. 

Now assume r is written in this canonical form; then ai ~ I and we have 
pq ~ I choices for the element 02 = ''"('^i) following ai in the cycle, pq — 2 
choice for the next element ~ ''"(02), and so on until the final element Op of 
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the first cycle. Then starting element Op+i of the second cycle has to be the 
minimum element of L \ {oi, 02, . . . , a^}, but we have {p — l)q — 1 choices for 
cip+2 = Tio,p+i)' an iterative argument proves the assertion. □ 

Lemma A. 3. The number of distinct self-conjugate permutations on a set of I 
elements is given by {I — 1)!! ;= (/ — 1) • (/ — 3) • . . . • 1. 

Proof. A self-conjugate permutation r on a set L of I = 2m elements is the 
product of TO disjoint 2-cycles, and the the result follows from Lemma [A. 2 1 □ 



Combining Lemma lA. II and IA.3[ we immediately get the following estimate: 
(2m- 1)!! 



2to 



s$ N^ig.n) (2to- 1)!!, m = m^,i^,{g,n), 



where the upper bound comes from assuming that no two clovers can be trans- 
formed one into the other by a rotation, and the lower bound comes from con- 
sidering all clovers as part of the same equivalence class. 

In order to create a trivalent fat graph from a clover, we replace the vertex 
with a full binary tree with I = 2m leaves; equivalently, we identify the leaves 
of the tree according to the same "gluing pattern" that identifies half-edges in 
the clover. 

More precisely, let Gq — (£; ctq, (Xi, 0-2) be a clover, with L = {!,..., 2to} 
and (To = (12 . . . 2m) as above. Let L' be set of leaves of a chosen binary tree T 
and / : L' ^ L a bijection. Now r := f~^ o o / is a fixed-point free involution 
on L': by identifying leaves of T according to r, we get a trivalent fatgraph G, 
which we say is obtained by plugging T into Gq (by means of /). 

Given a permutation (/)' on L' , the map /' = / o ^' is a bijection and we 
have: 

r' = /' ^ o ai o /' = ^ o {f~^ o (Ji o /) o 0' = 0' ^ o t o 0', 

which is an involution on L' conjugate to r. Conversely, if a'^ ~ <p~^ o ai o (/) is 
conjugate to cri, then f'~(f>of: L' — i is again a bijection, hence: 

/-I o a'l o / = (/-I o 0-1) o ai o (0 o /) = f'-' o ai o /', 

which is the involution defining the attachment map of T to Gq by means of 
/'. Since any two involutions arc conjugate, we can fix the map / once and for 
all binary trees with the same number of leaves, and only let the involution cti 
(i.e., the clover Go) vary. 

Therefore N2{g, n) = ^2(5, n) ■ F(m„iin((7, n)), where Y(l) is the count of full 
binary trees with / leaves. The number Y{1) is given by the (I — l)-th Catalan 
number: 

Hence from (|A.2|) we get: 

N^{g,n) ^ N2{g,n) N+{g,n), 
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where: 



N:^{g,n) 



1 (4m -2)! 

2^ ' (2TO-2)!!(2m)!' 



(2i-"'2)!i(L)r («> 

m := mnunig,n). 



A. 3 Generation from permutations 

As in the previous section, represent a fatgraph G in combinatorial form as 
(L; CTo, (Ti, cr2) where L = {1, . . . , 2m}. Here we count the number of trivalent 
fatgraphs that are generated by naively constructing a fatgraph from its com- 
binatorial definition. 

If G is trivalent, then cfq is a product of disjoint 3-cycles; by Lemma IA.2I 
the number of such ctq is: 

C(3,fc) = - - 2) • (; -4)(/ - 5) • ... -2 • 1, l = 2m = 'Sk (17) 

For each chosen ctoj each choice of a self-conjugate permutation ai gives rise 
to a trivalent (g, n)-fatgraph; by Lemma |A.3I there are exactly (2m — 1)!! such 
choices. Therefore, we have: 

Nsig.n) = (2m- 1)!! ■ C(3,2m/3) = (2m- 1)!! • (2m- l)(2m-2)- 

■(2TO-4)(2m-5)-...-2-l, ''^^^ 

where m — TOmax(3, n). 



B Pseudo-code notation 

Blocks of code arc marked by indentation (rather than delimited by specific 
keywords) . 

The 'def keyword is used to mark the beginning of a function definition. 

The notation 'for x in S" is used to loop over all the items x in a set or 
sequence 5*; sometimes the notation 'for x in a, is used instead. The form 

'for I, X in enumerate{Sy is used for keeping track of the iteration number when 
looping over the elements of S: as x runs over the items in S, i orderly takes 
the values 0, 1, . . . , up to \S\ — 1. 

B.l Basic types 

Numbers and basic data structures (arrays, lists, sets; sec below) are considered 
basic types, together with the logical constants True and False, and the special 
value None. 



B.2 Objects 

The word "object" is used to denote a kind of aggregate type: an object is a 
tuple '(fli, a2, . . ., Cat)', where each of the slots can be independently assigned 
a value; the values assigned to different a^'s need not be of the same type. We 
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write X.tti <— b to mean that the slot of object X is assigned the value 
b. Unless otherwise noted, object slots are mutable, i.e., they can be assigned 
different values over the course of time. 

An objects' class is the tuple '(ai, . . ., ajv)' of slots names that defines the 
object; the actual tuple of values is called an instanceV^ 

B.3 Arrays, lists, sets 

A few types of basic data structures are used in the code: arrays, lists and sets. 
They are distinguished only for clarity, and we make no assumption that these 
are primitive: for instance, each of these data structures could be implemented 
on top of the "list" type defined here. 

An "array" is a fixed-size collection of elements of the same type; the number 
and type of elements stored in an array will be stated when the array is first 
created. Items in an array can be accessed by position: if a is an array, then its 
fc-th element will be accessed as a[k]. Array elements can be mutated; we write 
a[k] <— b to mean that object b is stored into the fc-th place of array a. 

A "list" is a variable-size collection of objects. Two features distinguish lists 
from arrays: (1 ) lists can grow and shrink in size, and (2) lists can store items 
of different types. If Z is a list with n elements, the notation I. appendix) will 
be used to mean that x should be added as {n + l)-th item in list /. Again, the 
square bracket notation I [k] is used to denote the value stored in the fc-th place 
in I, and l[k] <— a; means that the fc-th slot of I is mutated to the value x. The 
operator "-I-" stands for concatenation when applied to lists. 

A "set" is a mutable unordered collection of objects of the same type. The 
only relevant difference with sets in the mathematical sense of the word is that 
set variables are mutable: if s is a set, then s.add{x) will be used to specify that 
s should be mutated into the set s U {x}. No duplicates arc admitted: it x G s 
and X = y, then s.add{y) does not alter s in any way. 

The word "sequence" will be used to denote any one of the above three. 
When S* is a sequence, we define size{S) as the number of elements in S; if S is 
a list or array object, valid indices into S range from to size{S) — l. 

B.3.1 List comprehensions 

A special syntax is used to form a list when its items can be gotten by applying 
a function or operation to the elements of another sequence. 

The notation 'L <r- [f{x) for x in S*]' makes L into the list formed by eval- 
uating function / on each element in S, analogously to the usual notation 
{f{x) : X & S} for sets. 

As an extension, the expression 'L <— [f{x) for x in S if P{x)y makes L into 
the list of values of / over the set S" of elements of S for which the predicate 
P(x) is true: S' = {x : x e S A P(x)}. 

-'^'^ Readers familiar with object-oriented programming will note that this is an over-simplified 
version of the usual object-oriented definition of objects and classes; this originates in the fact 
that the concrete implementation of the algorithms was done in object-oriented Python, but 
object-orientation is by no means essential to the implementation. 



58 



B.4 Operators 



The "%" operator is used to take the remainder of integer division: for integers 
k and n > 0, the expression (k % n) evaluates to the smallest non-negative 
residue of k mod n. 

The "+" operator normally denotes addition when applied to numbers, and 
eoncatenation when applied to lists. 

Any other operator keeps it usual mathematical meaning. 
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